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PREFACE 


The final report of Project RSC-3458, "Measurement of Soil Mois- 
ture Trends with Airborne Scatterometers" is divided into three vol- 
umes. The first volume deals primarily with the work completed by 
Dr. Sidney Theis relating multi spectral (visible through microwave) 
information to soil moisture trends in bare and vegetated fields. 
The second volume deals primarily with the work of Dr. Wesley 
Rosenthal in relating the same multi spectral data sets to agricultural 
crop classification and biomass estimation. The third volume by Ms. 
Cheryl Jones, details field work, aircraft schedules, data processing 
and calibrations, and the final data sets. 



ABSTRACT 


Monitoring soil moisture changes by remote means requires timely 
and repetitive sensor coverage. This dictates that the sensor must 
have all-weather capabilities. Microwave sensors have this capabil- 
ity. Until recently, microwave remote sensing of soil moisture has 
been limited to determining the best single sensor with the maximum 
sensitivity to soil moisture. Yisible/near-infrared sensors have been 
developed that estimate biomass and discriminate between crops. Vege- 
tation over a particular area does not change as rapidly as surface 
soil moisture. Therefore the absence of all-weather capabilities is 
not a severe limitation when visible/near-infrared biomass estimates 
are combined with microwave sensors to estimate soil moisture. 

Multi frequency sensor data from NASA's C-130 aircraft were used 
1) to determine which of the all-weather microwave sensors demon- 
strated the highest correlation to surface soil moisture over optimal 
bare soil conditions, and 2) to develop and test techniques which use 
visible/infrared sensors to compensate for the vegetation effect in 
this sensor's response to soil moisture. Soil moisture and aircraft 
data were collected at two agricultural areas: Guyrnon, Oklahoma in 
1978 and Dalhart, Texas in 1980. 

The L-band passive microwave radiometer was found to be the most 
suitable single sensor system to estimate soil moisture over bare 
fields. In comparison to other active and passive microwave sensors 
the L-band radiometer 1) was influenced least by ranges in surface 
roughness, 2) demonstrated the most sensitivity to soil moisture 
differences in terms of the range of return from the full range of 
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soil moisture, and 3) was less sensitive to errors in measurement in 
relation to the range of sensor response to soil moisture than the 
active microwave systems. L-band emissivity related more strongly to 
soil moisture when moisture was expressed as percent of field 
capacity. The perpendicular vegetation index (PVI) as determined from 
the visible/infrared sensors was found useful as a measure of the 
vegetation effect on the L-band radiometer response to soil moisture. 
A linear equation was developed to estimate percent field capacity as 
a function of L-band emissivity and the vegetation index. The predic- 
tion algorithm improves the estimation of moisture significantly over 
predictions from L-band emissivity alone. 



INTRODUCTION 


Estimation of soil moisture over large areas is useful in meteor- 
ology, hydrology, and agriculture. In meteorology, atmospheric models 
require information about the energy flux at the earth's surface. The 
two types of energy exchanged at the surface are sensible heat and 
latent heat. Sensible heat absorbed and released by the soil is for 
the most part, a small component of the surface energy balance. When 
latent heat is considered, some measurement of soil wetness is needed 
In order to relate actual evaporation to potential evaporation rates 
(Priestley and Taylor, 1972}. In hydrology, real-time estimates of 
Soil moisture condition at the beginning of a storm event would 
improve the ability to estimate runoff and provide flood warnings that 
would save both lives and property. In agriculture, a time history of 
surface soil moisture can be used to determine the severity and areal 
extent of drought conditions as well as inputs into soil moisture 
profile models needed to estimate agricultural yields. This 
information is useful to monitor the food supply for the growing 
population of the world. 

Until the advent of the computer and space age technology there 
was no reasonable way to handle the volume of data required to map 
soil moisture estimates over an area as large as the United States, 
much less on a global basis. Conventional soil moisture measurement 
techniques are very time-consuming and therefore measurements are not 
widely or regularly obtained over most of the United States. The 
spatial variations of soil moisture make it difficult to extrapolate 
conventional point measurements to represent an integrated value over 



a large area. These restrictions have led to research in developing 
remote sensing techniques to estimate soil moisture. The remote sens- 
ing approach to soil moisture estimation requires a new scale of 
thought. Sensing devices, particularly when satellite mounted, inte- 
grate many of the micro-scale variations that have been noticeable in 
point measurements. However, these remotely sensed measurements may 
actually give the operational user more meaningful information than 
point measurements, when used with a new set of models which use 
remote sensor inputs (Blanchard, 1979). 

To monitor surface soil moisture by remote means, timely and 
repetitive sensor coverage is required. Soil moisture changes rapidly 
in the near-surface layers. After an irrigation or rain, the soil 
surface dries out within a few days in clear sunny weather (Taylor and 
Ashcroft, 1972). This dictates that the sensor must have all-weather 
capabilities. Visible and infrared sensors do not meet the all- 
weather criterion because the short wavelengths in this region of the 
spectrum are attenuated and many times completely masked by cloud 
cover. Microwave sensors with wavelengths greater than 4 cm are not 
hindered by cloud cover or lack of sunlight. Also, longer wavelengths 
penetrate into the soil medium to a certain extent (Lundien, 1966). 
This provides microwave sensors with two potential advantages over 
visible and infrared sensors: 1) all-weather capability, and 2) the 
ability to gather information about the subsurface as well as surface 
moisture conditions. 

At visible and infrared frequencies, the sensor response contains 
information only about a very thin layer at the air-soil interface 
(Ulaby, 1974). However, this does not mean that visible and infrared 
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sensors are not valuable tools for large area soil moisture monitor- 
ing. Techniques using visible/infrared data have been developed that 
discriminate between crop types and estimate biomass, leaf area index, 
and percent cover (Richardson and Wiegand, 1977 and Rouse et_ al . , 
1973). Jackson et jil_. (1981b) suggested the use of vegetative biomass 
estimates from visible/near-infrared measurements in combination with 
microwave sensors to determine soil moisture. Vegetation over a 
particular area does not change as rapidly as surface soil moisture. 
Only three or four scenes in the visible and near-infrared wavelength 
region of the spectrum may be necessary to classify the vegetation 
during a season. Therefore the absence of all-weather capabilities is 
not a severe limitation when visible/infrared biomass estimates are 
combined with microwave sensors estimates of soil moisture. 

Three basic terrain parameters affect microwave sensors: 1) soil 
moisture, 2) roughness, and 3) vegetation. The frequency and incident 
angle determine the scale of roughness to which microwave sensors are 
sensitive. The influence of vegetation is also dependent upon fre- 
quency and incident angle. Until recently, research in microwave 
remote sensing of soil moisture has been limited to determining the 
best single sensor which maximizes sensitivity to soil moisture while 
minimizing the effects of roughness and vegetation. However, the 
optimum frequency and incident angle may change for each agricultural 
region with different tillage practices and types of crops that are 
grown. 

Active microwave investigations at the University of Kansas have 
identified the 4.25 GHz active microwave with incident angles near 10 
degrees as the one sensor that can measure soil moisture over all 
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terrain conditions (Ulaby and Batlivala, 1976; Ulaby jat aK, 1978). 
However, data presented by Ulaby et_ aU (1977) showed a 3 dB shift in 
the sensor response and a reduction in sensitivity caused by rough- 
ness. Active microwave has also shown promise in crop classification 
when used at large incident angles where the effects of roughness and 
vegetation dominate (Bush and Ulaby, 1978; Schwarz and Caspall, 1968; 
El lermeler et al . , 1967). 

Passive microwave sensors are also affected by the terrain para- 
meters. Numerous investigations, such as Newton (1977), Schmugge 
(1978, 1980a, 1980b, and 1980c), and McFarland (1976) have demon- 
strated the sensitivity of passive microwave sensors to soil mois- 
ture. Lee (1974) and Newton (1977) demonstrated that surface rough- 
ness influences this sensitivity to soil moisture. The effect of 
vegetation on passive microwave response was identified in studies 
such as Lee (1974) and Ki rdiashev ert jj1_. (1979). 

Very few investigations have combined the responses from various 
sensors to estimate soil moisture. Using L-band passive microwave 
system, Newton (1977) used the response from two polarizations at one 
incident angle to infer the roughness contribution on the microwave 
response at a different incident angle. Investigations which combine 
microwave sensors with sensors from other regions of the spectrum have 
not been performed because concurrent data have not been available. 

The attention of this study is focused on investigating techniques 
which infer soil moisture over agricultural lands from a combination 
of aircraft-mounted sensing devices. The sensors included visible, 
infrared, and microwave. Multi -sensor data over a wide range of 
frequencies and incident angles are available only from NASA's C-130 
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aircraft. Two of the most comprehensive soil moisture experiments 
using this aircraft were performed jointly by Texas A&M University and 
the University of California at Santa Barbara at sites in Guymon, 
Oklahoma and Dalhart, Texas. These data were used in this study 

1) to determine which of the all-weather microwave sensors 
demonstrates the highest correlation to surface soil moisture over 
optimal bare soil conditions. 

2) to develop and test techniques which use visible/infrared sensors 
to compensate for the vegetative effect in the microwave sensors' 
response to soil moisture. 

It is hypothesized that the final technique would have the potential 
for estimating soil moisture over a variety of agricultural crops with 
only remote sensing inputs. These techniques could then be expanded 
to larger areas with a satellite mounted sensor package. 

PHYSICS OF MICROWAVE INTERACTIONS 

There is some conceptual difficulty in relating to the microwave 
portion of the spectrum. Visible and near-infrared sensors are 

extensions of our eyes. Thermal sensors can be identified with 
walking barefoot over hot pavement or cool, lush grass. However, we 
cannot put our hand on a surface and feel the parameters which 

influence microwave response. The emission, transmission, scattering, 
and absorption of microwave radiation are governed by the same 
physical laws that govern radiation in the visible and infrared 

portions of the spectrum. Microwave energy penetrates through the 

atmosphere and clouds with minimal attenuation, giving the microwave 
sensors all-weather capabilities. Since microwave energy has longer 
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wavelengths, they are affected by parameters below the surface as well 
as those on the surface. This affords the opportunity to gather 
information about the soil volume itself. 

The ability of microwave sensors to detect soil moisture is based 
upon the difference in dielectric constant of dry soil and of water. 
It is common practice to express the dielectric constant, k, in terms 
of the dielectric constant of free space, k 0 : 

K p = k q K (1) 

where Kp is a dimensionless complex function expressed by 

< r = K r ' - j Kp " (2) 

where Kp' is the relative permittivity, < r " is related to conduc- 

tivity, and j denotes the complex portion of the function. The magni- 
tudes of both <p' and x r " are much greater for water than for dry 

soils and air space at microwave frequencies. In a dry soil, the 
value of Kp' is typically between 3 and 5, whereas the value of 
Kp' is about 80 for water (Schinugge et a]_. , 1974). 

Investigations such as those by Cihlar and Ulaby (1974), Hoekstra 
and Delaney (1974) and Newton (1977) have presented the relationship 
between permittivity (kp' and Kp") and moisture content for 

several different soil-texture types and frequencies. An example of 
the texture effects at 1.4 GHz is presented in Fig. 1 (Newton, 1977). 
Both the sand and clay soils have approximately the same value of 
Kp' at 0% moisture. The function of Kp' for each soil has two 

regions separated at a "transition" value of water content. For sand 
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FIG. 1 The relationship between relative permittivity and volumetric water content at 1.4 GHz for 
(a) sand and (b) Miller clay. From Newton (1977). 





the transition value is approximately 3 % moisture content in contrast 
to 20% for the clay. Newton (1977) attributes the change of < r ' 
below the transition zone with the increase of the tightly bound or 
hygroscopic water content. As more water is added, the more loosely 
bound pore water causes a larger effect (slope). This shift agrees 
with the established principle that with increasing clay content there 
is an increasing percentage of tightly bound unavailable water (Brady, 
1978). Newton (1977) concludes that the dielectric properties are 
directly related to availability of soil water to vegetation indepen- 
dent of soil type. It is important to note that < r is a function of 
frequency. 

Microwave sensors are of two basic types, active and passive. 
Active microwave sensors (radar and scatterometers) illuminate the 
target with microwave energy and measure the amount of energy returned 
to the antenna. Passive microwave sensors (radi ometers) measure the 
amount of microwave energy that the target emits. 

The terrain-sensor interaction processes in microwave remote 
sensing are governed by the target geometry and dielectric proper- 
ties. These properties depend upon their relation to the sensor 
parameters (Cihl ar and Ulaby, 1974). Sensor parameters include 
frequency, incident angle, and polarization as well as whether the 
sensor is active or passive. The ability of microwave sensors to 
detect soil moisture is obtained from the direct relationship between 
microwave response and the dielectric properties of the soil (Hoekstra 
and Delaney, 1974). As the soil moisture changes, its dielectric 
properties also change. This causes a corresponding change in the 
microwave response. 
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Target geometry refers to the slope, aspect, degree of roughness, 
and type and extent of vegetative cover (biomass). The effects of 
target parameters are all integrated into one microwave response. 
Consequently, it is difficult to relate physical measurements of 
target parameters to the effects they have on microwave responses. 
This is further complicated because these effects are also a function 
of the sensor parameters (Rouse et al_. , 1969). 


Active Microwave 

Radar return is that portion of the transmitted energy which 
returns to the radar receiver. The radar cross section, o, or alter- 
nately the radar cross section per unit area, (scattering coefficient) 
0 °, depends upon the properties of the target. Thi sparameter contains 
all the information about the illuminated terrain that a radar is 
capable of sensing. The scattering cross section, o m , for a single 
point scatterer is related to the radar return, by inverting the radar 
equation and solving for a m (Rouse et al . , 1969): 


o 

m 



( 4 ">V 
^ **• 


( 3 ) 


Where: 

o m = scattering cross section for a single point scatterer 

0 " 2 ) 

P r = received signal power (watts) 

P-t = transmitted signal power (watts) 

D m = slant range distance to the mth scatterer (m) 

\ = wavelength of the transmitted signal (m) 
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G m = antenna gain (di mens i onl ess ) 

The average scattering coefficient per unit area for many scatterers 
is a dimensionless real number whose magnitude is a function of the 
terrain parameters. Since a 0 is dimensionless it can be analyzed in 
dB and expressed as a function of terrain parameters: 

o° = f ( x, e, <)>, P, k, r) (4) 


where 

e = angle of incidence relative to nadir 
<j> = aspect angle 

P = polarization of the incident wave 
k - complex dielectric constant 
F = index of surface roughness. 

An analytical expression for o° is not readily obtainable without 
simplifying assumptions (Rouse et al_. , 1969). 

In the following discussion, a perfectly smooth surface at the 
air-soil interface is assumed to facilitate the understanding of the 
relationship of a 0 to soil moisture. Under these conditions the 
Fresnel reflection coefficient, r, can be defined. This coefficient 
relates the magnitude and phase of the reflected electric and magnetic 
fields to those of the incident fields of a plane wave which strikes 
the interface. The parameter of interest is the power reflection 
coefficient, R: 



(5) 


R H 


F(p-u) 2 9 + qj~ 1 

[(p+y) 2 + q 2 ] 


( 6 ) 
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( 7 ) 


_ [(>V ~ p)2+ ( *V - q) 2 ] 

v ' [(< + P) 2 + (<’ - q?] 

where H and V represent horizontal and vertical polarization and p = 
cose. The parameters p and q are given by: 


D - 1 

P — 1 ’ 

KV 

. 2 . \ 2 , H 1 1 /2 

+ P - 1) + K r J 

+ <K r ' + M 2 - 1)) 1/2 

(8) 

q = 72~ 

<[<v' 

■ 2 i\2i H 1 1/2 

+ P - 1) + Kp ] 

- U r " + M 2 - 1)} 1/2 

(9) 


The sensitivity of R to the complex dielectric constant is not readily 
apparent from these equations. Using approximate values of < r from 
Fig. lb and equation (6), R^ for a range of look angles is presented 
in Fig. 2. Small incident angles show the most sensitivity to changes 
in soil moisture. The term "look angle" is synonymous with incident 
angle. 

The Fresnel approach addresses only smooth or specular return 
from the surface with no subsurface contribution or depolarization. 
The energy is reflected as from a mirror, with the angle of incidence 
equal to the angle of reflection. Spectral reflection is seldom pres- 
ent in natural- terrain except for areas with smooth water surfaces. 
At the other extreme, isotropic scatterering is when the scattering 
coefficient is independent of the angle of incidence. The "degree of 
roughness" of the terrain, a relative quantity, dictates the extent of 
scattering which forms the reradiation pattern. Examples of specular 
and isotropic scattering along with an intermediate case are presented 
in Fig. 3. 
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FIG. 2 Relationship between power reflection coefficient (R^) and 
look angle for 0% and 30% soil moisture. 
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3 Examples of spectral and isotropic scattering along with an 
intermediate case. Incident angle is defined by e. 
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A large number of scattering theories have been developed by 
Barrick (1968), Kodis (1966), Fung (1966), Beckmann (1965), Semenov 
(1966), Beckmann and Spizzichino (1963), Leader and Dalton (1972), 
Rouse (1972) and Blanchard (1977). The authors have modeled the 
effects of surface roughness, subsurface contribution, and 
depolarization. These modeling efforts have had varying degrees of 
success for particular terrain roughnesses but no one theory works 
universally over all scales and distributions of roughnesses and 
subsurface conditions. 

Vegetation effects further complicate the dependency of c° on 
soil moisture. Vegetation scattering models have been developed by 
Peak (1959a), Du (1969), and Waite and Cook (1974) with varying 
degrees of success. The physical phenomena of radar scattering is 
very complex and not well understood. The scope of this study does 
not allow in-depth investigations of radar reflectivity theory. Such 
a study must be done on small, controlled surface areas with accurate 
measurements of roughness, vegetation, and soil dielectrics. However, 
it is important to recognize the impact of scattering from surface 
roughness and vegetation on the scattering coefficient in the attempt 
to relate the scattering coefficient to soil moisture. 

Passive Microwave 

Since the microwave emission of a soil originates in the soil 
volume and propagates outwardly, it is reasonable to examine the emis- 
sion process within the volume separately from the modifying effects 
of surface roughness and vegetative cover. The fundamental concepts 
of microwave radiation are presented here. These concepts are des- 
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cri bed in more detail by Paris (1969 and 1971), Marion (1975), Newton 
(1977), Hess (1959), and Lintz and Simonett (1976). 

The specific intensity of radiant energy emitted by a blackbody 
is expressed by Planck's radiation law: 


"bb =1!l ? [exp (h v/CT )-l ] -1 
X 


( 10 ) 


where 

Bbb a specific intensity of radiant energy emitted in a 

particular frequency range that passes through an element 
of area in an interval of time and surrounded by a solid 
angle (watts nr 2 steradiarrl Hz~l) 

C = Boltzmann's constant (1.38 x 10~ 23 JK~1) 

T = temperature of the blackbody (K) 

A ■ wavelength (in) 

h = Planck's constant (6.623 x lO'" 3 ^ Js) 
v = frequency (Hz) 

At microwave frequencies and surface temperatures typical of 
Earth, the Rayleigh-Jeans approximation of Planck's equation for 
thermal radiation from a blackbody applies. A Maclaurin expansion of 
the bracketed quantity in equation (10) is: 


exp (h v/CT )-l = hv/CT + h 2 v 2 /2C 2 /T 2 + ... 


( 11 ) 


Equation (10) reduces to the Raylei gh-Jeans approximation 


R « 2CT 
B bb “7? 

A 


( 12 ) 
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when the higher order terms in equation (11) are neglected as very 
small . 

The emission from a blackbody is completely random; therefore, 
the polarized specific intensity of emission (Bp) can be expressed 
as 


B 


P 



(13) 


since one-half the total intensity of emission would exist in each 
plane of polarization. Material objects on the earth's surface are 
not perfect radiators, so they are sometimes referred to as gray- 
bodies. The intensity of graybody radiation may be thought of as the 
intensity radiated by an equivalent blackbody multiplied by an effici- 
ency factor. It is assumed that the form of the emission is the same 
as that of a blackbody and that the polarized intensity, Ip, is 
given by 

Ip=cB p (14) 


where e is the polarized ernissivity of the gray 
The polarized brightness temperature, Tg-j-, 


body. 

is defined by 



which reduces to 

T bt - Cl 


(15) 


(16) 
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where T is the temperature of the emitting graybody. Simply stated, 
all objects with temperatures above absolute zero emit electromagnetic 
energy due to the thermal agitation of charged particles. The radia- 
tion intensity is linearly related to the temperature and physical 
properties of the object. 

Various sources contribute to the radiation received by a passive 
microwave system. The most important contribution for this study is 
the contribution from the earth's surface (gTsurf)* Other major 
sources include the atmosphere between the antenna and the surface, 
T/\tm> and downward radiation from both sun and the atmosphere that 
is reflected from the surface to* the radiometer, Tsky* For a 
constant height, frequency, and angle 

T BT = T ATM + T ( ( 1 - e ) r sKY + e T SURF^ 

where 

T( 3 t = brightness temperature, as seen by the radiometer (K) 
t = atmospheric transmission (dimensionless) 

1-e = reflectance of the surface (dimensionless) 
e = einissivity of the surface (dimensionless) 

TsuRF = thermoinetric temperature of the surface layer (K) 

Most studies, including this one, neglect all atmospheric effects 
(Tatm, T^ky=0, t«1). This simplifies equation (17) to 

T BT = e T SURF‘ 
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This - assumption is valid for the majority of atmospheric conditions, 
but as shown by Paris (1971), it loses validity in the presence of 
relatively large particles such as rain drops in the atmosphere. 

When applied to a soil this simple approach to microwave emission 
assumes a semi -i nf inite soil medium with a smooth surface that is in 
thermodynamic equilibrium and exhibits blackbody (isotropic) radia- 
tion. Peake (1959b) demonstrated that e is equal to absorptivity. 
This is the basis for describing the emission from a soil volume, e 0 , 
in terms of the power reflection coefficient: 

e 0 = 1 - % ( 19 ) 

where p denotes polarization and Rp can be calculated from equation 
(6) or equation (7). In Peake's analysis uniform soil moisture and 
temperature distributions with depth were assumed. Richerson (1971) 
and Casey (1972) used an effective Fresnel coefficient from horizontal 
homogeneous layers to provide for non-uniform moisture profiles. 
Strogryn (1970) developed a rigorous solution for non-uniform 
temperature and moisture profiles using Maxwell's equations,, Using an 
approach similar to Strogryn, Tsang et_al_. (1975) formulated emission 
from a stratified medium. Numerical results were illustrated and 
compared with the closed form analytical solutions. 

Radiative transfer techniques (Newton, 1977; Burke and Paris, 
1975) provided a simpler alternative to the general formulation of 
Stogryn (1970). The soil volume was modeled as a horizontally strati- 
fied media with each layer consisting of non-scattering homogeneous 
soil. Soil moisture and temperature were assumed constant in each 
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layer,, This approach lends itself to physical interpretation and is 
not computationally difficult. 

In nature, the soil surface is seldom perfectly smooth and with- 
out vegetation. The remainder of this section examines the influence 
of surface roughness and vegetation upon the passive microwave res- 
ponse to soil moisture. 

Surface roughness modifies the microwave emission from a soil 
volume. Choudhury et_ _aL (1979) modeled roughness by modifying 
equation (19) with a roughness parameter, h: 

e = Gq + Ae (20) 

where 


U - R p (l - e' h ) (21) 

Using both 1.4 GHz and 10.6 GHz radiometers, Lee (1974) and Newton 
(1977) reported measurements of bare, smooth, and rough soil surfaces 
that demonstrated that increased surface roughness decreases the sen- 
sitivity of the emission to changes in soil moisture, and that for dry 
soils the emission is approximately independent of surface roughness. 
Using these same frequencies, B 1 inn and Quade (1974) obtained measure- 
ments for both raked and smooth sand. From their results, the 10.6 
GHz frequency was affected much more by that specific scale of 
roughness than the 1.4 Ghz frequency. This indicates that the 
"sensed" roughness is a function of wavelength. For a smooth field 
Newton (1977) reported a 100 K change in Tgj in going from wet to 
dry soils. This sensitivity was reduced to about 60 K for extremely 
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rough fields that had been deep plowed to produce large clods. This 
shows good agreement with the 80 K range of Tbt for the moderately 
rough fields encountered in Phoenix (Schmugge, 1978). 

The effects of vegetation on the surface brightness temperature 
have been modeled by Sibley (1973), Basharinov and Shutko (1978), and 
Kirdiashev et_ aj_. (1979). Experimental results from Lee (1974) sup- 
port Sibley's model which predicts that the vegetation is essentially 
an attenuator at low vegetation density and a predominant emitter at 
high vegetation density. As an attenuator, the vegetation lowers the 
apparent temperature measurement; as a predominant emitter, it raises 
the apparent temperature. Lee (1974) used 1.4 GHz and 10.6 GHz 
frequency radiometers and reported that for 10.6 GHz, thick vegetation 
had completely masked the soil contribution regardless of surface 
roughness. In contrast, for the 1.4 GHz frequency, the vegetation 
generally exhibits an attenuating effect. Therefore, the effect of 
vegetation on apparent temperature is dependent on frequency, with 
lower frequencies showing more potential to detect soil moisture 
through the vegetation. 

REVIEW OF PREVIOUS WORK 
Active Microwave 

Early ground-based measurement programs at Ohio State University 
(Taylor, 1959 and Cosgriff et_ aj_., 1960) and by the Corps of Engineer 
(Lundien, 1966) investigated radar responses and penetration into a 
variety of natural and man-made surfaces. Roughness, vegetation, and 
soil moisture were found to be the general classes of parameters which 
influenced radar response. Due to the wel 1 -developed photographic 
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interpretation techniques, the greatest excitement in the remote 
sensing community during the 1960s and early 1970s was created by 
radar image analysis. Numerous authors presented qualitative image 
interpretations of the effects of vegetation (Ellermeier et aL , 1967; 
Moore and Simonett, 1967; Morain and Simonett, 1966; Schwarz and 
Caspall, 1968; and Berger, 1970) and suggested radar as a crop 
classifier. However, the launching of Landsat redirected vegetative 
studies away from radar and toward visible/infrared techniques. 

MacDonald and Waite (1971) renewed interest in radar's potential 
for soil moisture detection. Their study was the first to qualita- 
tively discern soil moisture variations in radar imagery at relatively 
high frequencies and large incident angles. Previous experiments with 
similar imaging radars had been conducted with soil moisture showing 
little influence upon radar return even for bare soils (Simonett, 
1970). However, MacDonald and Waite (1971) discerned the difference 
between wet swampy lands and the better drained, relatively drier 
soils of natural levees. Both these areas were beneath defoliated 
tree branch cover. Studying flooded low areas in Central Texas, 
Blanchard et _al_. (1979) described this phenomena as multiple scatter- 
ing between the trees and the smooth free water surface. 

The University of Kansas has since provided the bulk of research 
into the dependency of radar response on soil moisture. A description 
of the truck-mounted Microwave Active Spectrometer (MAS) system is 
provided in Ulaby et_aj_. (1978). The research using this system iden- 
tifies C-band radar (4.25 GHz) at incident angles near 10 degrees as 
the optimum frequency which can estimate soil moisture over all ter- 
rain conditions with a minimum of influence from roughness (Ulaby and 
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Batliyala, 1976; Ulaby etjiK, 1978). However, sensitivity (approxi- 
mately 3 dB) to roughness (Ulaby et_ _al_. 1977) and vegetation type 
(Ulaby _et^ _al_. » 1979) is demonstrated at this frequency and incident 
angle. They found that as soil moisture varies from 0 to 30%, o° 

changes by approximately 10 dB. Jackson et _al_. (1981a) reported 
agreement with the University of Kansas with aircraft mounted radars 
by illustrating that C-band at 15 degrees incident angle exhibited the 
best linear correlation with soil moisture. However, the ).4 points 
presented exhibited a sensitivity of over 20 dB per 30% moisture 

range. 

Blanchard (1978) presented results from ai rcraft-mounted radar 
experiments (13.3 and 1.6 GHz). The influence of roughness and 

vegetation were minimized at angles between 10 and 20 degrees. At the 
10 degree incident angle, the moisture sensitivity increased from 4 dB 
at 13.3 GHz to 9 dB at 1.6 GHz per 30% soil moisture. This closely 

agrees with the Kansas data for 4.25 GHz reported by Ulaby et_ a]_. 

(1977) and indicates no significant difference in sensitivity between 
4.25 and 1.6 GHz. However, the 1.6 GHz data from Blanchard (1978) 
does not agree with the low sensitivity shown by Ulaby at 1.1 GHz. 

Most radar measurements are accurate to only 3 dB. As mentioned 
before, roughness has at least a 3 dB effect at the optimum incident 
angles. The magnitude of these variations along with vegetation 
effects can significantly mask the effectiveness of radar to quantify 
soil moisture when applied to a variety of terrains. 
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Passive Microwave 


Many investigators have experimentally shown the sensitivity of 
truck and aircraft-mounted passive microwave radiometers to soil 
moisture variations. Poe et al_. (1971) observed that the emissivity 
of a smooth bare field varied from 0.5 for very wet soil to greater 
than 0.9 for dry soil. Blanchard (1972) and Blinn and Quade (1972) 
reported decreasing Tbj with increasing moisture content. The rates 
ranged from 1.5 to 4.8 K per one percent change in moisture content 
depending on frequency and soil type. This range is supported by 
Basharinov _et a]_. (1974)', who reported a slope of approximately 3-4 K 
per one percent moisture change for the 1-3 cm wavelength range. A 
comparison of aircraft observation of Tbj over agricultural fields 
around Phoenix in 1973 (Schmugge et al . , 1974) and 1975 (Schmugge, 
1976; Choudhury et_ aL , 1979) demonstrated that these results are 
repeatable. Recent controlled truck-mounted experiments by Newton 
(1977) and Wang et_ aj_. (1980) agree closely with the previous 
findings. 

Schmugge (1980a) developed a relationship between soil texture 
and field capacity 

FC = 0.30 - 0.0023 x SAND + 0.005 x CLAY (22) 

where: 

FC = field capacity (for use with volumetric soil moisture) 

SAND = percent sand • 

CLAY = percent clay. 
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He determined that passive microwave response should be related to 
percent field capacity instead of soil moisture by volume. Percent 
field capacity may, in some applications, be more important than 
absolute water content. 

The sensitivity of radiometers to soil moisture has also been 
demonstrated from space. McFarland (1976) showed a definite 
relationship between the Skylab 21 cm brightness temperatures and an 
antecedent precipitation index (API). The soil moisture index (API) 
was calculated from rain gage data within the 110 km footprint. 
Eagleman and Lin (1976) also related Skylab 21 cm data to soil 
moisture variations as determined by water budget models. The 1.55 cm 
Electrically Scanning Microwave Radiometer (ESMR) response was 
correlated to moisture indices (McFarland and Blanchard, 1977; Theis 
1979; Schmugge et_ ah, 1977). This shorter wavelength radiometer 
demonstrated the potential to detect soil moisture during periods of 
the year when the agricultural fields were essentially bare. 

Visible/Infrared Biomass Models 

Reflected solar radiation is the primary source of energy for 
visible/infrared sensors. Different materials possess different spec- 
tral reflective properties unique to that material. Due to Fresnel 
reflectance at air/water interfaces within leaves, radiation at the 
near- and middle-infrared frequency is strongly reflected (Gates, 
1980). This reflectance increases with increasing biomass. This 
phenomenon has been the basis for classification of vegetation and 
estimating biomass. Two biomass estimation models are presented here. 
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Richardson and Wiegand (1977) developed the perpendicular vegeta- 
tion index (PVI). This graphical approach presented in Fig. 4 plots 
two Landsat bands against each other; MSS5(0.6-0.7 pm) and 
MSS7 (0. 8-1.1 pm). As soil reflectance properties change, the two 
bands change together to form the soil line. PVI was calculated using 
the Pythagorean theorem; 

PVI = \/(Rgg5 - Rp5) 2 + (Rgg7 - Rp7) 2 (23) 


where Rp is the reflectance values of an individual band and Rgg is 
the corresponding soil reflectance for the Landsat bands. PVI is the 
perpendicular distance from the soil background line to the plotted 
radiance values. 

The Transformed Vegetation Index (TVI) has been used as a range- 
land biomass estimator (Rouse et a]_. , 1973); 


TVI 


(MSS7 - MSS 5 
: MS'S 7 + MS'S 5 


+ U. 5 


(24) 


The ratioing under the radical is a standard practice to remove the 
scene to scene variations of sun angle and atmospheric optical density 
difference. The 0.5 was added to prevent the term under the square 
root from becoming negative. PVI was used in this study because it 
was developed for agricultural crop lands and may be used as both a 
crop classifier and biomass estimator. 

Sensor Combinations 

Very few investigations have combined responses from various 
sensors to estimate soil moisture. Newton (1977) used the difference 
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FIG. 4 Diagram illustrating the principle of the perpendicular 

vegetation index (PVI) model. A perpendicular from candidate 
plant coordinates (Rp5, Rp 7) intersects the soil background 
line at coordinates (Rg5, Rg7). A PV 1=0 indicates soil, and 
a PVI>0 indicates vegetation. From Richardson and Wiegand 
(1977). 
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in the horizontally and vertically polarized L-band radiometer at 35 
degree incident angle as a parameter to infer roughness. This 
parameter was used in combination with the L-band emissivity at 20 
degree incident angle as an estimate of soil moisture over various 
roughnesses. 

Jackson et al. (1981b) developed a model to predict emissivity of 
a scene from surface temperature, soil moisture, and vegetative bio- 
mass inputs. Encouraging results were obtained for wet biomass values 
below 1300 g/m 2 . Beyond this, the sensitivity is probably too low for 
the model to be of use. These results support the need for biomass 
estimations (possibly obtained from visible/infrared sensors) as an 
indicator of vegetation effect on bare soil emission. Jackson could 
not investigate this because visible/infrared data were not taken 
simultaneously with the microwave data. The data collected for this 
study were used to develop sensor combination techniques to estimate 
soil moisture by remote means. 

DATA COLLECTION AND PROCESSING 

Multisensor data over a wide range of frequencies are available 
only from NASA's C-130 aircraft. The C-130 collected all sensor data 
used in this study. Aircraft and soil moisture data were collected 
near Dalhart, Texas in August 1980 and near Guymon, Oklahoma in 
August, 1978. These areas were chosen because each had relatively 
homogeneous soils, reasonably low relief, a variety of crops, and the 
availability of irrigation. The Texas and Oklahoma Panhandle areas 
were initially surveyed in one mile square sections with boundaries of 
the sections oriented North-South and East-West. Aircraft flight 
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lines were also orientated in these directions to make it possible for 
the pilot to fly one half mile over the fields that were sampled and 
thus allow adequate sampling by the sensors over each field. 
Agricultural practices usually dictate one crop per quarter section 
which is also the area irrigated by one center pivot sprinkler 
system. For this study a field was defined as one half of a quarter 

section with the long axis of the field parallel to the aircraft 

flight line. The half quarter section provided an area large enough 

for sufficient coverage with aircraft line sensors and small enough 

for an adequate ground sampling network. 

Experiment Site Descriptions 

Guymon 

Sampled fields were located under four flight lines approximately 
20 km southwest of Guymon, Oklahoma (Fig. 5). The soil type was 

generally a silty clay (35% clay, 35% silt, 'and 30% sand) with many 

areas having caliche (CaC0 3 ) near the surface. The crop types were 

alfalfa, and milo, with some fallow fields. The milo fields were 

chosen with some fields having rows perpendicul ar and some parallel to 
the flight lines. The field conditions are summarized in Table 1. 

Dalhart 

The sampled fields were chosen under two flight lines approxi- 
mately 20 km northwest of Dalhart, Texas (Fig. 6). The soil type of 
the surface 30 cm was a sandy loam (75% sand, 10% silt, and 15% 
clay). The 22 fields sampled were almost equally divided into 

irrigated and non-i rri gated fields. The irrigated fields included 
millet and mature corn. Pasture, wheat stubble, and disked wheat 
stubble comprised the non-i rrigated fields. 
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FIG. 5b Legend for Guymon, Oklahoma field maps. 
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Table 1. Surface Conditions of Guymon Fields 


Condition 

Fields Numbers 



Summer Fallow (bare) 

2, 10, 17, 21, 26 

Sprinkler Irrigated August 1 

14 

Sprinkler Irrigated August 14 

6 

Listed August 17 

2X 

Circular Irrigated Alfalfa 

4, 13, 22, 27 

Milo 

Circular Irrigated 

23 cm tall, drilled 

7, 15 

90 cm tall perpendicular rows 

20, 25 

Furrow Irrigated 

90 cm tall parallel rows 

8, IX, la , 2a 

90 cm tall perpendicular rows 

18, 24 


Table 2. Surface Conditions of 

Dalhart Fields 


Condition 

Field Numbers 


Bare 

Wheat stubble 

Disked stubble 

Mulched stubble (with weeds) 

Half disked, half burned stubble 

• 17,18 
19,20,21,22 
13,14 
15,16 


Corn 

1,2,7,8,9,10,11,12 


Pasture 

5,6 


Mi 1 let 

3,4 
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[G. 6 Dalhart field maps, a) Area map of Dalhart showing relative 
location of each field map. 
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FIG. 6e Locations of sample fields at Dal hart 




The summer of 1980 vas extremely dry In the Dal hart area. 
However, It was fortunate that this area received a one Inch tain 
after the first sampling date. This provided a range of soil moisture 
conditions. Those fields with vegetation were reasonably uniform In 
crop cover. All corn fields were In mature stages. Only fields 1 and 
2 showed areas of stress especially where some of the sprinkler 
nozzles were not functioning properly. Field conditions are 
summarized in Table 2. 


Soil Moisture Data 

Soil moisture samples were collected at eight points per 32 
hectare field (Fig. 7) except in the Guymon fields 6 and 14 which were 
intensively sampled (37 points). Dalhart sample depths were 0-2 cm, 
2-5 cm, 5-15 cm, 0-15 cm, 15-30 cm, 30-45 cm (Fig. 8). The 5-15 cm 
interval was divided into two samples at Guymon; 5-9 cm and 9-15 cm. 
Manpower limitations did not allow each field to be sampled during the 

To correlate the ground data to aircraft data, the soil moisture 
values in fields not sampled on flight days had to be extrapolated. 
To do this, graphs were generated plotting soil moisture versus time. 
By incorporating precipitation events and Irrigation schedules in the 
graphs, it was possible to produce dry-down curves for each field and 
point. This method also enabled the estimation of soil moisture for 
points where values were missing due to pesticide spraying, rain, 
sample spills, etc. Where the curve for each depth intersected a 
flight time, the value was recorded as the new gravimetric soil 
moisture. These values plus the bulk density results were used to 
calculate the volumetric soil moistures by point. Texture analyses 
were performed for each field in Guymon and for selected fields in 


Dalhart. 
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FIG. 7 Sampling pattern for fields at Guymon and Dal hart 
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FIG. 8 Soil moisture sampling depths at Dal hart and Guymon. The 
15-30 and 30-45 cm core samples were also taken in addition 
to the above. The 5-15 cm sample was divided into two 
samples at Guymon; 5-9 cm and 9-15 cm. 
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Ground tracks of the line sensors were used to determine which 
sample points fell within the area covered by the sensors. These 
points were arithmetically averaged by field and termed soil moisture 
under the line sensors or just soil moisture. 

The soil sample processing lab was set up with conventional 
microwave ovens. As the samples came in from the field, they were 
weighed with seals and lids on the cups. This weight (measured to a 
tenth of a gram) was recorded as the wet weight. Lids and seals were 
then discarded and soil samples dried in the ovens. At regular time 
intervals two test samples per oven were removed and weighed and 
returned to the ovens. The weights were plotted on a graph to track 
the dry-down. When the weight loss of the test samples leveled out to 
less than a 0.2 gram difference from the previous weighing, they were 
considered dry. All samples were then removed from the oven and their 
dry weights recorded. As the preprocessing was completed for a set, 
the samples were stored for possible future reference. 

A programmable calculator was used in the lab to compute gravi- 
metric soil moistures as the final weights were completed. The pro- 
gram automatically subtracted an average lid, seal and cup weight from 
the gross values. Daily soil moisture charts were kept on each field, 
point and depth. This system facilitated quality control of the 
data. If a value was questionable, the sample was located, reweighed 
and its soil moisture recalculated immediately. 

Aircraft Data 

Aircraft sensor data were similar for both experiments. The 
Guymon data were collected on August 2, 5, 8, 11, 14, and 17, 1978. 
The Dalhart data were collected on August 14, 16 (two flights), and 18 
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(late afternoon), 1980. The sensor package consisted of 1) seven 
scatterometer frequencies and polarizations, 2) three passive 
microwave radiometer frequencies and polarizations, 3) Barnes PRT-5 
thermal data, 4) aerial photography, and 5) visible/infrared scanner 
data. Soil moisture flights were made at an altitude of 500 m. 
Scatterometer frequencies and polarizations included: 

1) 13.3 GHz VV 

♦ 

2) 4.75 GHz HH 

3) 4.75 GHz HV 

4) 1.6 GHz HH 

5) 1.6 GHz HV 

6) 0.4 GHz HH 

7) 0.4 GHz HV 

where VV = vertically polarized transmitted and received 

HH = horizontally polarized transmitted and received 
HV = horizontally polarized transmitted and vertically received. 
Eight look angles from nadir were processed: 5, 10, 15, 20, 25, 35, 
40, and 45 degrees. Near-nadir (3 degree) passive microwave data were 
collected at 1) L-band (1.6 GHz) horizontal polarization, 2) C-band 
(4.75 GHz) vertical, and 3) C-band horizontal polarization. In 
Dal hart, the eight channels of the NS001 (simulated thematic mapper) 
bands included: 

1) Channel 1: 0.45-0.52 pm 

2) Channel 2: 0.52-0.60 pm 

3) Channel 3: 0.63-0.69 pm 

4) Channel 4: 0.76-0.90 pm 

5) Channel 5: 1.00-1.30 pm 

6) Channel 6: 1.55-1.75 pm 
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7) Channel 7: 2.08-2.35 pin 

8) Channel 8: 10.4-12.5 pin 

In Guymon, five channels from the modular multi spectral scanner (M 2 S) 


were processed: 

1) Channel 4 0.548-0.583 pin 

2) Channel 7 0.662-0.701 pm 

3) Channel 8 0.703-0.747 M m 

4) Channel 9 0.770-0.863 pin 

5) Channel 11 8.000-12.00 pin 

The M 2 S channels 7, 9, and 11 correspond closely to the NS001 channels 
3, 4, and 8, respectively. 


Data Processing 

Processing procedures for the individual sensors and ground data 
are given in Appendix A. The procedures included calibration, quality 
control, data omission, and field average determination. Data omis- 
sion usually occurred with excessive roll and drift of the aircraft. 
Once field averages had been calculated for each sensor and their cor- 
responding soil moisture, the ground and aircraft data sets were 
merged. Further analysis was done on the field averaged data set. 

ANALYSIS AND DISCUSSION 
Approach 

Case studies were first performed to investigate the microwave 
sensor responses to soil moisture and roughness. Previous studies 
have related microwave sensor responses to surface soil moisture at 
near incident angles. Near incident angle scatterometer (active) and 
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radiometer (passive) microwave data for all available frequencies and 
polarizations were plotted and regressed against soil moisture (0-2 
cm) over essentially bare soils. Regression line slopes and 
correlation coefficients were used as criteria for evaluating the 
sensitivity and relative value of each sensor configuration. From 
this analysis a single microwave sensor, the L-band radiometer, was 
chosen as the sensor that exhibited the highest potential to quantify 
soil moisture. 

The analysis was then expanded to vegetated fields for this 
sensor. Linear relations were determined for the L-band response to 
soil moisture for each crop type. These individual relations were 
used to obtain predicted values of soil moisture. The correlation 

between measured and predicted soil moisture over all fields 
represented the optimum capability of the L-band passive microwave 
sensor to quantify soil moisture when crops can be accurately 
classified. 

In general, crop type is related to biomass; however, there is a 
range of biomass within each crop type. A visible/infrared biomass 
estimator and crop classifier (P V I ) was directly combined with the 
L-band passive microwave data to quantify soil moisture independent of 
crop type. The direct combination technique was compared to the 
classification techniques and microwave alone. 

Case Studies 

Two case studies were performed on the sensor data from Dal hart 
in order to investigate each microwave sensor's response to particular 
changes in 1) soil moisture and 2) roughness. Fig. 9 illustrates a 
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field being actively irrigated at the time of the data collection. 
Soil moisture measurements were not available but the field is notice- 
ably wetter (darker) on the left side where it had been irrigated 
recently. The field had uniform roughness on a scale similar to near- 
by sample fields. Line plots of each sensor are presented at the same 
horizontal scale as the photograph in Figs. 10-12. Passive microwave 
sensor (Fig. 10) exhibited approximately a 25 K sensitivity to the 
difference in soil moisture while the surface thermal temperature 
range was approximately 3.5 K. This small thermal response is due to 
the time of day the data were taken. The flight time was mid-morning 
but this change in thermal temperature indicates that there was a 
change in moisture between the two portions of the flight line. The 
scatterometers at 10 degree incident angle (Fig. 11) exhibited con- 
siderable signal noise and a maximum sensitivity of only approximately 
5 dB after the noise is smoothed. It is interesting to observe that 
all the scatterometers exhibit a stronger dependence on soil moisture 
at 40 degree incident angle (Fig. 12). This is contrary to results 
from most previous studies. 

The second case study involved field 15 at Dalhart. A photo is 
not presented but the left half portion of this field had been 
recently disked. The right half was burned wheat stubble (relatively 
smooth). The field had uniform soil moisture (15 percent by volume). 
A definite roughness influence is apparent in the passive microwave 
responses (Fig. 13) that does not readily appear in the 10 degree 
scatterometer data (Fig. 14). The roughness effect is again apparent 
in all the 40 degree scatterometer data (Figure 15) with the exception 
of returns from the 1.6 GHz HV, 0.4 GHz HH, and 0.4 GHz HV scattero- 
meters. 
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FIG. 10 Line plots (sensor response vs. time) of field in Fig. 9 for 
L- and C-band radiometers and PRT-5. 
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FIG. lib Line plots (sensor response vs. time) of field in Fig. 9 for 
1.6 GHz HH, 1.6 GHz HV, 0.4 GHz HH, 0.4 GHz HV scatterometers 
at 10 degree incident angle. 
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FIG. 12a Line plots (sensor response vs. time) of field in Fig. 9 for 
13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz HV scatterometers at 
40 degree incident angle. 
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FIG. 13 Line plots (sensor response vs. time) of case study field 
(rough and smooth) for L- and C-band radiometers and PRT-5. 
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SIGMA (2.54 cm = 5-dB ) 
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ROUGH SMOOTH 



RELATIVE TIME ( 2.54 cm = 2 sec ) 

FIG. 14a Line plots (sensor response vs. time) of case study field 

(rough and smooth) for 13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz 
HV scatterometer at 10 degree incident angle. 
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FIG. 14b Line plots (sensor response vs. time) of case study field 
(rough and smooth) for 1.6 GHz HH, 1.6 GHz HV, 0.4 GHz HH, 
and 0.4 GHz HV, scatterometer at 10 degree incident angle. 
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RELATIVE TIME ( 2.54 cm « 2 sec ) 


FIG. 15a Line plots (sensor response vs. time) of case study field 

(rough and smooth) for 13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz 
HV scatterometer at 40 degree incident angle. 
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FIG. 15b Line plots (sensor response vs. time) of case study field 
(rough and smooth) for 1.6 GHz HH, 1.6 GHz HV, 0.4 GHz HH, 
and 0.4 GHz HV scatterometer at 40 degree incident angle. 
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Microwave Sensor Responses to Soil Moisture Over Bare Fields 

This analysis is presented in two sections (scatterometers and 
radiometers) to avoid confusion between the active and passive sys- 
tems. For the same set of fields, plots of sensor responses versus 
volumetric soil moisture (0-2 cm) are presented for each sensor. The 
fields selected for these illustrations were all bare with a range of 
roughness effect due to tillage and surface weathering. There was no 
quantitative measure of roughness, but it was observed to range from 
smooth to moderately rough as defined by Newton (1977). 

Scatterometers 

As mentioned before, previous investigations indicated that the 
scatterometer 1 s highest potential for detection of soil moisture is 
around the 10 degree incident angle. Plots of each scatterometer's 
response at the 10 degree incident angle versus soil moisture in the 
surface 2 centimeters are presented (Fig. 16-19) to illustrate 
differences in the data sets. The scatterometer returns at low soil 
moisture shifted approximately 0 dB, 5 dB, 8 dS and 7 dB for the 13.3, 
4.75, 1.6 and 0.4 GHz systems, respectively. These shifts illustrate 
that the comparison of the two data sets must be approached with 
caution. Further analysis was conducted for each scatterometer within 
each data set to minimize the influence of lack of calibration. 

Each scatterometer' s response versus volumetric soil moisture 
(0-2 cm) was analyzed for 10, 15, 20, and 40 degree incident angles. 
The 1.6 GHz HH scatterometer (Dal hart) is illustrated in Fig. 20 as a 
representative sample of the scatterometer responses. The remainder 
of the illustrations are in Appendix B. The absence of a direct 
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FIG. 16 Scatterplot of 13.3 GHz a 0 vs. volumetric soil moisture (0-2 
cm) for Guymon and Dalhart bare fields. 
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FIG. 17 Scatterplot of 4.75 GHz o° vs. volumetric (0-2 cm) for Guymon and 
Dal hart bare fields. 
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LEGENO: SITE a □ a OALHART x x x GUTHON 

FIG. 18 Scatterplot of 1.6 GHz a 0 vs. volumetric soil moisture (0-2 cm) for Guymon 
and Dalhart bare fields. 
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LEGEND: SITE * □ a DRLHRRT X x x GUTHON 

FIG. 19 Scatterplot of 0.4 GHz a 0 vs. volumetric soil moisture (0-2 cm) for Guymon 
and Dal hart bare fields. 
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and soil moisture for all bare fields and each scatterometer (10, 15, 
20 and 40 degree incident angles) are compiled in Table 3. 


Table 3. Scatterometer R 2 Values for o° versus 
Volumetric Soil Moisture 



R 2 Values 



Guymon 

OaU 

lart 


Scatterometer 

10° 

15° 

o 

O 

CVJ 

o 

o 

10° 

15° 

o 

O 

CVJ 

o 

13.3 

VV 

mm 

0.47 

0.65 

0.10 

0.20 

0.29 

0.21 


4.75 

HH 


0.35 

0.46 

0.60 

0.13 

0.14 

0.12 


4.75 

HV 

HH 

0.36 

0.49 

0.51 

0.45 

0.49 

0.45 

0.36 

1.6 

HH 

0.48 

0.58 

0.47 

0.35 

0.27 

0.46 

0.43 

0.28 

1.6 

HV 

0.43 

0.59 

0.57 

0.62 

0.02 

0.28 

0.36 

0.17 

0.4 

HH 

0.02 

0.11 

0.19 

0.22 

0.08 

0.38 

0.19 

0.07 

0.4 

HV 

0.02 

0.00 

0.00 

0.00 

0.29 

0.49 

0.22 

0.17 

The 

largest 

R 2 va 

lue (0 

.65) obtained 

for either 

data 

set was 


Guymon's 13.3 GHz VV (20 degree); however, the largest average R 2 


value (0.52) was obtained for the 1.6 GHz HH (15 degree). These 
values will be used in the next section to compare the capability of 
the scatterometers to detect soil moisture to that of passive 
microwave radiometers over the same set of fields. 


Radi ometers 

Radiometers provide intensity measurements in terms of brightness 
temperature (K). By assuming the temperature of the emitting layer to 
be that sensed by the PRT-5, emissivity can be calculated using 
equation (18) as 

e = T BT/ T PRT-5 ( 25 ) 
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For this study the ernissivity was used as the measure of the passive 
microwave response. The data sets were initially analyzed together 

(Figs. 22-24) to determine if there were calibration problems. No 

significant calibration problems were evident except in the L-band 
horizontal ernissivity (Fig. 24) where the dry ernissivity values were 
significantly higher at Dalhart; however, the linear R 2 value for this 
combined data set was 0.71. The radiometer data were further analyzed 
by comparisons to volumetric soil moisture and percent field capacity 

for each data set (Figs. 25-27). Both the C-band horizontal and 

vertical data showed considerable scatter for Dalhart (Figs. 25 and 

26). The scatter is related to the differences in the roughness of 
the field groups. Fields D17&18 appeared "smoother" to the 
radiometers than fields D15&16. Recall that these two fields also 
exhibited corresponding differences in relative roughness in the 40 
degree incident angle 1.6 GHz HH scatterometer (Fig. 20). The same 
roughness effect is apparent in the L-band radiometer data but to a 
lesser degree (Fig. 27). The "sensed" relative roughness did not 
influence the longer L-band wavelengths as much as it did the C-band 
wavelengths. 

Percent field capacity (PFC) is defined as 

PFC = VSM/FC x 100 (26) 

where 

VSM = measured volumetric soil moisture 

FC = Field capacity obtained from texture analysis and equation 

( 22 ) 


67 


'</)«/>— in jn< 02330 


paul o 

OF POOR QUALITY 


DALHART AND GUYMON BARE FIELDS 
C-BAND VER. EMISSIVITY VS SOIL MOISTURE 


\x °x 

x \°r 

\&jp p 

*x*bV x 

v ilr 3 x 

x x 

xCfa \ 


R 2 =0.53 


I 0.85- 
V 


x ’ 
x * X 
□ 


LECEN'Ds H 


SOIL MOI ST'J^E tO- 2 CM) 


xxx OSLHSr.T a a o G'JTMON 


FIG. 22 Scatterplot and regression line for C-band (vertical 

polarization) emissivity vs. volumetric soil moisture for 
Dal hart and Guymon bare fields analyzed together. 
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FIG. 23 Scatterplot and regression line for C-band (horizontal 

polarization) emissivi.ty vs. volumetric soil moisture for 
Dal hart and Guymon bare fields analyzed together. 
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FIG. 24 Scatterplot and regression line for L-band (horizontal 

polarization) emi'ssivity vs. volumetric soil moisture for 
Dalhart and Guymon bare fields analyzed together. 
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FIG. 26 C-band (horizontal polarization) response to volumetric soil 
moisture and percent field capacity for Guymon and Dalhart 
bare fields analyzed separately. 
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FIG. 27 L-band (horizontal polarization) response to volumetric soil 
moisture and percent field capacity for Guymon and Dalhart 
bare fields analyzed separately. 
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The difference in soil texture between the Guy in on and Dal hart test 
sites is better identified in the L-band radiometer measurements since 
the effect of roughness, which could be confused with texture effects 
is reduced at this frequency. Guymon had a much higher clay content 
than Dalhart. The slopes of the regression lines are different at 
Guyinon and Dalhart when e is related to volumetric soil moisture 
(-0.0074 and -0.0011, respectively). In contrast, they are almost 
identical when e is related to percent field capacity. This supports 
the hypothesis that passive microwave responses should be analyzed as 
a function of percent field capacity. After consideration of the 
field capacity plots it is evident that the calibration problem with 
L-band was magnified in the percent field capacity plots because it 
was not offset by the texture effect. The difference in the intercept 
between data sets is 0.043. This correction was added to the Guymon 
data set for all further analyses. The correction was added to Guymon 
because the L-band antenna used in collecting that data set was 
deteriorating and replaced between the times of the Guymon and Dalhart 
experiments. The deterioration tended to reduce the passive microwave 
response at Guymon (Blanchard, 1981). 

The one inch rainfall which occurred at Dalhart provided an 
opportunity to relate the roughness "sensed" by L-band radiometers 
with that sensed by L-band scatterometers . This was not possible in 
Guymon because the irrigated bare fields (G6&14) provided only one 
field with a sufficient moisture range to establish the relative 
roughness effects. The o° intercept values at 0 percent moisture, 
o° 1 , were obtained for each of the Dalhart bare field groups from the 
1.6 GHz HH (40 degree) scatterometer (Table 4). 
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Table 4. o 0 ' Values for Each Field Group in Dalhart 


Field 

1 

c o 

D13&14 

-30.8 

D15&16 

-28.7 

D17&18 

-36.2 

D19&20 

-32.6 

D20&21 

-.31.5 


This "roughness" was considered to influence both the intercept and 
the slope of the e vs percent field capacity relationship. The linear 
equation: 

e = A + B * PFC (27) 


was modified by 

A = A1 + A2 * o o ‘ ' (28) 


B = B1 + B2 * o 0 ‘ (29) 

e = A1 + (B1 8 PRF) + (B2 * PFC * a Q ' ) + (A2 * a ') (30) 

The resulting equation when applied to bare fields in the Dalhart data 
set is: 

e = 1.212 + (0.0016 5*PFC ) + (0.000145*PFC* a °' ) + (.00631 x a 0 ') (31) 

with R 2 value = 0.93 and F value = 158.5. The significant increase in 
R 2 values indicates that scatter in the L-band radiometer is indeed 
due to difference in relative roughness as sensed by the 1.6 HH (40 
degree) scatterometer. From the case study and comparison of R 2 
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values, the L-band radiometer shows the strongest relationship to soil 
moisture over the same set of bare fields than any other microwave 
sensor. Therefore, L-band radiometer data were used in all further 
analyses. 

L-band Emissivity Response to Soil Moisture over all Crop Types- 

No significant difference in e response to moisture could be dis- 
cerned when certain "near bare" fields were included in the classifi- 
cation of bare fields. These included Guymon's 23 cm tali mi 1 o fields 
(7 and 15) and Da 1 ha rt ' s pasture (5 and 6) and millet (3 and 4). This 

indicates that for L-band, attenuation of the soil emission is insig- 

nificant for these fields. For this experiment the crop types were 
classified and generally can be rank-ordered from lower to higher 
bi omass: 

1) bare soil 

2) alfalfa 

3) milo 

4) corn 

Four fields of alfalfa were sampled at Guymon. Fields 4 and 13 
and fields 22 and 27 were pairs of adjoining fields. During the 18 

day experiment the alfalfa was irrigated and matured. Fields 22 and 

27 were harvested before the last flight day. These two fields were 
on a steep hillside; consequently, the effective incident angle was 
increased from 3 degrees to greater than 10 degrees. An increase in 
effective incident angle in the horizontal polarization causes a 
decrease in emissivity (Newton, 1977). The decrease was observed in 
Fig. 28a where the values of e for fields G22&27 were generally 
shifted downward from the values for fields G4&13. A realistic 
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correction of 0.02 (Newton, 1981) was added to the emissivity for 
fields G22&27 with a resulting increase in the significance of the 
linear relationship between e and PFC (Fig. 28b) (R 2 values increased 
from 0.49 to 0.67). All mi 1 o fields considered in this study were at 
the Guymon site while all corn fields were at the Dalhart site. The 
linear analysis of the mi 1 o fields (Fig. 29) showed a strong relation- 
ship between e and PFC (R 2 value = 0.69). The same analysis for corn 
fields (Fig. 30) yielded a much weaker relationship (R 2 value = 0.08) 
due to: 

1) Corn masked or attenuated the emissivity from the soil, thus 
destroying the response to soil moisture differences. 

2) The difference of biomass depending on the field group varies 
markedly within the corn classification. 

The influence of different biomass within the corn fields is 
noted when PVI values are compared to regression lines for individual 
field groups. Fields D7&3 and D11&12 had PVI values greater than 
4.5. These fields were combined into field D>4.5. Regression lines 
for the three subclasses of corn are presented in Fig. 31. The 
relatively flat slope for the combined data for high biomass corn 
fields (D>4.5) indicates that as PVI values approach 4.5 the measured 
emissivity is virtually insensitive to the soil moisture under the 
corn. To 'illustrate the need to compensate for vegetation, emissivity 
was plotted versus percent field capacity for all fields combined with 
the resulting R 2 values of 0.09 (Fig. 32). 

Most previous studies have modeled emissivity as the dependent 
variable while using soil moisture and biomass as independent 
variables (Jackson at al . , 1981b; Choudhury, 1981). Choudhury used 


the model 
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FIG. 29 Scatterplot and regression line for L-band emissivity vs. 
percent field capacity for Milo fields. 
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FIG. 30 Scatterplot and regression line for L-band emissivity vs. 
percent field capacity for Corn fields. 
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FIG. 31 Scatterplot and regression lines with corresponding PVI 

values for L-band emissivity vs. percent field capacity for 
individual fields of corn. 
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e,„ * E s e" T + (1 - e $ ) e y (l - e" 1 ) e" T + Ey (l - e' T ) 


( 32 ) 


where t = vegetation index 

e m = measured emissivity 

e s = emissivity of the soil (function of soil moisture) 
e v = emissivity of vegetation 

The first term on the right is the contribution from the soil surface 
as attenuated by the vegetation. The second term is the contribution 
from the vegetation as reflected by the soil surface (l-e s ) and 
again attenuated by the vegetation. The third term is the direct 

contribution from the vegetation. By rearranging terms: 


'm 


— 2 T / *"* T 

= e y - e y e + e s ( e " e - 


e" T + e e“ 2-r ) 


( 33 ) 


and solving for e s , the resulting equation is: 


, 0 ~ 2 t 

. c in ~ e v + E e v 

Cs ’ e' T - e e~ T + e e' 2T 
V V 


( 34 ) 


There is always an error in the measurement of e m . was 

replaced with e m + e er (e er is the error in measurement) 
yielding 


e s = 


e + e 
m er 


+ e e 
v 


e ~ T - Ey e e y e - 2t 


( 35 ) 


In a similar fashion, the error term can be applied to the emissivity 
of the soil e s resulting in the following equation 
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e m = e y - e v e” 2T + (e s + e er )(e“ T - e v e" T + e y e" 2x ) (36) 

The soil moisture information is contained in 55. Choudhury 

reported values of e v to be 0.97. Using the value e e p 3 0.01, the 

sensitivity of es and e r to different values of t was investigated 
in Fig. 33. It was noted that as t increases from 0 to 1.5 the error 

in eg increases from 0.01 to 0.18 while the error of q n decreases 

from 0.01 to 0.0006. This indicates that these equations can predict 
the emissivity measured at the radiometer well, but cannot estimate 
the es or emissivity of the soil accurately when the vegetation 
becomes dense. 

Coefficients of the linear regression equations for the calcula- 
tion of e from percent field capacity for each crop type are tabulated 
in Table 5. 

Table 5. Coefficients of linear prediction equations for 
estimation of emissivity as a function of percent 
field capacity. 


Crop Type 

Intercept 

Slope 

Bare 

0.9801 

-.00253 

Alfalfa 

0.9638 

-.00125 

Mi 1 0 

0.9869 

-.00136 

Corn 

0.9657 

-.000121 


These coefficients were used to predict the value of emissivity. The 
calculated values of emissivity resulting from these equations are 
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alMAU PAGE 


Illustration of the sensitivity of Choudhury's model to 
errors in measurement vs. a vegetation index (t). The curve 
denoted by e s is the error in the predicted emissivity of 
the soil that would be incurred with a 0.01 error in the 
measured emissivity. The curve denoted by e m is the error 
in the predicted measured emissivity with a 0,01 error in 
e s* 




compared to the measured emissivity in Fig. 34. The \\ z value of 
measured versus predicted emissivity using this approach is 0.74. 

Soil Moisture Determination with Multiple Sensors 

The remainder of this study considers the soil moisture parameter 
as the dependent variable. Regression lines for percent field capac- 
ity versus emissivity for the four major crop types along with mean 
vegetation index values are shown in Fig. 35. This figure introduces 
two approaches to estimate soil moisture by using visible/infrared 
data with passive microwave sensors. The first approach uses 
visible/infrared sensors to classify the crop and crop Glass to esti- 
mate moisture. The second approach bypasses classification by using 
the PVI values {obtained from visible/infrared sensor) and microwave 
response in a direct relationship. 

Classification Technique 

The classification technique approach assumes that all fields 
within a particular crop type have equal biomass and thus the same 
vegetation effect on the passive microwave radiometer. This assump- 
tion was tested by investigating the effect of different biomass (as 
indicated by PVI) within crop types for corn and milo. Regression 
lines for each field group with the corresponding average PVI values 
are presented in Fig. 36. Significantly different regression lines 
within both crop types were observed. The sensitivity of the percent 
field capacity to changes in emissivity was directly related to PVI. 
That is, as PVI Increased, the sensitivity of PFC to e increased. 
This increased sensitivity dictates that errors in the measurement of 
e have greater effects on the accuracy of estimates of the percent 
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FIG. 34 Scatterplot and recession line for measured vs. predicted 
emissivity for all fields. 
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FIG. 35 Regression lines relating percent field capacity to 

emissivity and the corresponding PVI values for each crop 
type. 
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FIG. 36 Regression lines relating percent field capacity to emissivity and the PVI values 
for each field group in a) milo and b) corn. 



field capacity at high PVl values. This same relation Is the control- 
ling element in the sensitivity analysis of the Choudhury equation 
that was discussed in the previous section. 

Coefficients for the linear regression equations used for esti- 
mating percent field capacity as a function of measured emissivity are 
presented in Table 6. The regression lines defined by these values 
were illustrated in Fig. 35, 

Table 6: Coefficients of linear prediction equations for 

estimation of percent field capacity as a 
function of emissivity 


Crop Type 

Intercept 

Slope 

Bare 

291.86 

-291.97 

Alfalfa 

493.61 

-493,65 

Milo 

512.70 

-510.19 

Corn 

707.31 

-656.58 


Measured percent field capacity was plotted against the predicted 
values obtained from these equations (Fig. 37a). This Illustration 
points out the failure of using one set of linear coefficients per 
crop type (especially in corn) to predict soil moisture for an entire 
set of fields inside that crop type. The most significant R 2 value 
which can be expected using classification technique on this data set 
with corn included is 0.67. When the corn fields were not considered 
the R 2 value increased to 0.76 (Figure 37b). Thus, if the severe 
attenuation from corn can be avoided we may be able to estimate per- 
cent field capacity from measured emissivity reasonably well. 
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MEASURED VS. PREDICTED PERCENT FJELO CAP 
(NON-CORN FIELDS) 
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Direct Combination 


The estimation of percent field capacity with a combination of 
sensors required the development of an equation that defines percent 
field capacity as a function of both microwave emlsslvity and Vegeta- 
tion index data (PVI). This technique should be more flexible than 
the classification technique described earlier since It can accommo- 
date differences in biomass within each crop class. 

The regression lines for the different crop types and the corres- 
ponding values of PVI are presented in Fig. 38 with the corn sub- 
divided Into three field groups. The slope's of these regression lines 
were then plotted against the average PVI values In Fig. 39. The 

form of the relation appears to be a straight line function when 

fields D9&10 and D>4.5 are not considered. The linear form of the 
relationship between percent field capacity and emisslvity is: 

PFC = A - B * e (37) 

If biomass (PVI) is considered to linearly affect both A and B then: 

A = A1 + A2 * PVI (38) 

B = B1 + B2 * PVI (39) 

Substituting into equation (37) and expanding the form, the equation 
becomes: 

PFC = A1 + (A2 * PVI) - (B1 * e) - (B2 * e * PVI) (40) 

Applying this to all fields except those with average PVI value 

greater than 4.3 (D9&10 and D> 4 . 5 ) yields the equation: 

PFC = 279.53 + (51.20 * PVI) - (281.22 * e) - (48.41 * e * PVI) (41) 
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FIG. 38 Regression Tines relating percent field capacity to 

emissivity for all fields (corn divided into individual 
field groups). 
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It is significant to note the values of A1 and B1 In equation 
(41) are comparable to the values of A and B obtained for bare fields 
(291.86 and 291,97, respectively) (Table 6). This Indicates that the 
influence of vegetation was compensated for by the Increased values of 
PVI obtained over the non-bare fields. Comparisons of the capability 
of the L-band radiometer to estimate percent field capacity are 
presented in Fig. 40 when It is used? 

1) alone 

2) 'with a classification technique which excludes all corn 

3) in direct combination with PVI excluding fields with average 
values of PVI greater than 4.3. 

This illustrates that passive microwave alone cannot estimate soil 
moisture when several crop types are considered. The classification 
technique works reasonably well when the variation of biomass within a 
crop type are small. This technique also requires perfect classifica- 
tion. The direct combination technique shows the greatest potential 
to detect soil moisture because it: 

1) is able to discern difference in biomass within crop types 

2) will vary as the crop matures 

3) is not hindered by mi sclassifi cation 

4) provides a means to determine a threshold value of biomass 
beyond which the radiometer is unable to accurately estimate 
soil moisture. 
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FIG. 40 The ability of the L-band radiometer to predict percent field capacity when used 

a) alone, b) with classification technique (no corn), and c) with direct combination 
technique (excluding fields with PVI greater than 4.3). 


CONCLUSIONS AND RECOMMENDATIONS 


The principal objective of this research was twofold. The first 
objective was to select the most useful single sensor from the array 
of sensors available for the estimation of soil moisture. The second 
objective was the development of a linear combination of sensor data 
that would provide an improvement in estimation of soil moisture over 
estimates made with a single sensor. 

The L-band passive microwave radiometer was found to be the most 
suitable single system. Several subordinate conclusions led to the 
selection of this sensor. First the L-band system was found to be the 
most sensitive to soil moisture differences in terms of the range of 
return from the full range of soil moisture. Secondly, the error in 
measurement in relation to the range of sensor response was small in 

"r 

relation to the active microwave systems. Third, it became obvious 
that with the passive microwave systems it was possible to get a rela- 
tive calibration while with the active microwave systems this could 
not be done. Fourth, the L-band horizontal passive microwave res- 
ponded less to variations in roughness than the C-band passive micro- 
wave. In the preliminary analysis leading to these conclusions, it 
was shown that the L-band passive microwave measured emissivity rela- 
tionship to soil moisture was more significant when moisture was 
expressed as percent of field capacity. 

All previous experiments analyzed the microwave return, whether 
it was the passive microwave antenna temperature or the active micro- 
wave scattering coefficient, as a function of the moisture measure- 
ments. A more useful but more difficult modeling effort is required 


97 


to predict moisture from remote sensing Inputs. Hue to the nature of 
the Influence of surface roughness and vegetation on the transmission 
of microwave energy from the surface soil volume to the antenna, sig- 
nificant errors in estimated soil moisture will occur unless the sys- 
tems used can provide a reliable measure of surface roughness and 
vegetation attenuation. No practical satisfactory measure of rough- 
ness was found in these data; however, the perpendicular vegetation 
index was found to be useful as a measure of vegetation effects for 
agricultural crops other than dense corn. A linear equation was 
developed to predict percent field capacity as a function of L-band 
emissivity and the vegetation index. The prediction algorithm 
improved the estimation of moisture significantly over predictions 
from L-band emissivity alone. 

The results of this experiment infer that effective estimation of 
soil moisture will require as a minimum a combination of a dual polar- 
ized L-band radiometer (at some incident angle greater than 15 
degrees) and a visible/infrared imager with at least two bands that 
are suitable for developing the perpendicular vegetation index. Most 
importantly all of the sensors must have- a satisfactory means for 
calibration. An active microwave system could also be used for the 
roughness measurement if such a system is calibrated. Further experi- 
ments along these lines should be conducted with truck systems in 
order to maintain more control of the site and better calibration of 
instruments, while collecting and processing the data at lower cost. 


98 


REFERENCES 


1. Barrick, D. E. 19G8. Rough surface scattering band on the 
specular point theory. IEEE Trans, on Antennas and Prop. Ap-16 
(4 ) : 449-454 . 

2. Basharinov, A. Ye and H. M. Shutko. 1978. Uetermi nation of the 
moisture content of the earth's cover by superhigh frequency 
(microwave) radiometric methods: a review. Radio Engineering and 
Electronic Physics, 23:1-12, translation. 

3. Basharinov, A. Ye, L. F. Borodin and A. M. Shutko. 1974. 
Passive microwave sensing of moist soils. Preprints, 9th Int. 
Symp. Remote Sensing of Environ. April 15-19. Environmental 
Research of Michigan, Ann Arbor, MI. p. 3G3-367. 

4. Beckman, P. 1965. Scattering by composite rough surfaces, 
Proc. IEEE 53:1012-1015. 

5. Beckman, P. and A. Spizzichino. 1963. The scattering of 
electromagnetic waves from rough surfaces. Macmillan Co. New 
York. 

6. Berger, D» H, 1970* Texture as a discriminant of crops on radar 
imagery, IEEE Trans., Geosci. Elect. GE-8(4) : 344-340. 

7. Blanchard, A. 0. 1977. Volumetric effects in the depolarization 

of electromagnetic waves scattered from rough surfaces. Ph.D. 

Dissertation. Dept, of Electrical Engineering, Texas ASM Univer- 
sity, College Station, TX. 192 p. 

8. Blanchard, B. J. 1972. Measurements from aircraft to charac- 

terize watershed, 4th Annual Earth Resources Program Review, 
NASA, Manned Spacecraft Center, Houston, TX. Jan, 17-21, 

V (1 IB) : 1 -14. 

9. Blanchard, B. J. 1973. Measurement of soil moisture trends with 

airborne scatterometers. Progress Rep. 3458-2, Remote Sensing 

Center. Texas A&M University, College Station, TX. 102 p. 

10. Blanchard, B. J, 1979. Soil moisture. Proc. Fifth Annual 
William T. Pecora Memoml Symp. on Remote Sensing. American 
Water Resources Association, Sioux Falls, $D. June 10-15, p. 
331-332. 

11. Blanchard, B. J. 1980. Measurements of soil moisture trends 
with airborne scatterometers. Progress Rep. RSC 3458-4. Remote 
Sensing* Center. Texas A&M University, College Station, TX. 67 p. 

12. Blanchard, B. J. 1981. Personal Communications. 


99 



13. Blanchard, 13. J., 0, L. Nieber and A. J. Blanchard. 1979. 

Continuation of measurements of hydrologic soil cover complex 
with airborne scatteroineters. Final Rep. 3496. NSG-5156, Remote 

Sensing Center, Texas A&M University, College Station, TX. 43 p. 

14. Blinn, J, C. and J. G. Quade. 1972. Microwave properties of 
geological materials: studies of penetration depth and moisture 
effects. 4th Annual Earth Resources Program Review NASA Manned 
Spacecraft Center, Houston, TX. Jan. 17-21. 11:53-1 to 53-12. 

15. Brady, N. C. 1978. The nature and properties of soils. 8th 
edition McMillan Publ. Co. New York. 639 p. 

16. Burke, W. J. and J, F. Paris, 1975. A radiative transfer model 
for microwave emission from bare agricultural soils. NASA Tech. 
Memo. TMX-58166. NASA/JSC. Houston, TX. 

17. Bush, T. F. and F. T. Ulaby. 1978. An evaluation of radar as a 
crop classifier. Remote Sensing Environ. 7:15-36. 

18. Casey, K. F. 1972. Application of Hill's Functions to problems 
of propagation in stratified media. IEEE Trans, on Antennas and 
Propa. AP-20(3) : 368-374. 

19. Choudhury, B. J. 1981. Personal Communication. 

20. Choudhury, B. J. , T. J. Schmugge. R. W. Newton, and A. Chang. 
1979. Effect of surface roughness on the microwave emission from 
soils. J. Geophys. Res. 84:5699-5706. 

21. Cihlar, J. and F. T. Ulaby. 1974. Dielectric properties of 
soils as a function of moisture content. RSL Tech. Rep; 177-47. 
The University of Kansas Center for Research, Lawrence, Kansas. 

61 p. 

22. Clark, B. V. and R. W. Newton. 1979. JSME scatterometer data 

processing. Final Report 3337. NAS9-14875, Remote Sensing 

Center, Texas A&M University, College Station, TX. 110 p. 

23. Claassen, J. P., R. 0. Stroud, B. V. Clark, B. R. Jean and R. W. 
Newton. 1979. The system and hardware design of real-time fan 
beam scatterometer data processors. Final Report 3556, Remote 
Sensing Center, Texas A&M University, College Station, TX. 275 p. 

24. Cosgriff, R. L.y W. H. Peake and R. C. Taylor, 1960. Terrain 
scattering properties for sensor system design (terrain handbook 
II) Engineering Experiment Station. College of Engineering, The 
Ohio State University, Columbus. Bulletin Vol. 29, No. 3, 118 p. 

25. Dellwig, L. F. and R. K. Moore, 1966. The geological value of 
simultaneously produced like- and cross-polarized radar imagery, 
J. Geophys. Res. 71(14) : 3597 -3601 . 


100 


26 . 


Du, L. 1969. Scattering and absorption by large leaves at 

microwave frequencies. Tech. Rep. 2440-6 Electro. Science 
Laboratory. Ohio State University, Columbus, OH. 

27. Eagleman, J. R. and W. C. Lin. 1976, Remote sensing of soil 

moisture by a 21 cm passive radiometer, J. Geophys. Res. 81*. 
3660-3666. 

28. Ellermeier, R. 1)., D. S. Simonett and L. F. Dellwig. 1967, The 
use of multi -parameter radar imagery for the discrimination of 
terrain characteristics. IEEE , International Convention Record, 
Part 2, New York, p. 127-135. 

29. Fung, A. K. 1966. On depolarization waves backscattered from a 
rough surface. Planetary Space Sciences 14:563-568. 

30. Gates, D. M. 1980. Biophysical ecology. Springer-Verlag. New 

York. 611 p. 

31. Hess, S. L. 1959. Introduction to theoretcal meteorology. 

Holt, Rinehard and Winston. New York. 355 p. 

32. Hoekstra, P., and A. Delaney. 1974. Dielectric properties of 

soils at UHF and microwave frequencies. 0. Geophys. Res. 79: 
1699-1709. 

33. Jackson, T. J. , A. Chang, and T. J. Schmugge, 1981a. Aircraft 

active microwave measurements for estimating soil moisture, 
Photograin. Eng. and Remote Sensing 47:801-805. 

34. Jackson, T. J. , T. J. Schmugge and J. R. Wang. 1981b. Effects 

of vegetation on passive microwave estimates of soil moisture. 
Intern. Geosci. and Remote Sensing Symp., IEEE, Washington, D. 
C,, June 8-10, 1:375-387. 

35. Kirdiashev, K. P., A. A. Chukhlantsev and A. M. Shutko. 1979. 

Microwave radiation of the earth's surface in the presence of 
vegetation cover. Radiotecknika 1 Elecktronika 24: 256-264 

(translation) . 

36. Kodis, R. 1966. A note on the theory of scattering from an 

irregular surface. IEEE Trans, on Antennas and Propa. AP-14:77- 
82. 

37. Leader, J. C. and W. Dalton. 1972. Bidirectional scattering of 
electromagnetic waves from the volume of dielectric materials. 
J. Appl . Physics 43:3080. 

38. Lee S. L. 1974. Dual frequency microwave radiometer measure- 

ments of soil moisture for bare and vegetated rough surfaces. 
Tech. Rep. RSC-56. Remote Sensing Center, Texas A&M University, 
College Station, TX. 211 p. 


101 


39. Lint?,, J. and D. $. Sinionett. 1976. Remote sensing of 
environment. Addison-Wesley Publ. Co. Reading, Mass. 694 p. 

40. Lundien, J. R. 1966. Terrain analysis by electromagnetic means, 
Report 2, Radar responses to laboratory prepared soil samples, 
Technical Report No. 3-693, U. S. Army Corps of Engineers, U. S. 
Army Engineer Waterways Experiment Station, Vicksburg, Missis- 
sippi, September, 55 p. 

41. MacDonald, H. C. and W. P. Waite. 1971. Soil moisture detection 
With Imaging radars. Water Resources Research 7:100-110. 

42. McFarland, M. J. 1976. The correlation of Skylab L-band bright- 
ness temperatures with antecedent precipitation. Preprints 
Conf. on Hydrometeorology . Amer. Met. Soc. Ft. Worth, p. 60- 
65. 

43. McFarland, M. J. and B. J. Blanchard. 1977. Temporal correla- 
tions of antecedent precipitation with Nimbus 5 ESMR brightness 
temperatures. Preprints 2nd Conf. Hydrometeorology. Toronto, 
Ontario Canada. Amer. Meteor. Soc. p. 311-315. 

44. McFarland, M. J., C. Jones, S. W. Theis and W. D. Rosenthal. 
1981. Measurement of soil moisture with airborne scattero- 
meters. Progress Rep. RSC 345B-5. Remote Sensing Center. Texas 
A&M University. College Station, TX, 46 p. 

45. Marion, J. B. 1975. Classical electromagnetic radiation. 
Academic Press. New York. 479 p. 

46. Moore, R. K. and D. S. Simonett. 1967. Radar remote sensing in 
biology. Bioscience 17(6) : 384-390. 

47. Morain, S. A. and D. S. Simonett, 1966. Vegetation analysis with 
radar imagery. Proceedings, Fourth Symposium on Remote Sensing 
of the Environment, Report No. 4864-11-X, Willow Run Laboratories 
of the Institute of Science and Technology, The University of 
Michigan, Ann Arbor, April, p. 605-622. 

48. Newton, R. W. 1977. Microwave remote sensing and its applica- 
tions to soil moisture detection. Tech. Rep. RSC-81, Remote 
Sensing Center, Texas A&M University, College Station, TX. 500 p. 

49. Newton, R. W. 1981. Personal Communication. 

50. O'Neill, P. 1981. Personal Communication. 

51. Paris, J. F. 1969. Microwave radiometry and its application to 
marine meteorology and oceanography. Ref No. 69-11. Contract 
No. 2119(04). Dept, of Oceanography. Texas A&M University, 
College Station, TX. 210 p. 


102 



52. Paris, J. F. 1971. Transfer of thermal microwave in the 

atir.v.;ohere. NASA Grant NGR-44-001-098. Dept, of Meteorology 
Rpt . Texas A&M University, College Station, TX. 257 p. 

53. Peake, W. H. 1959a. Interaction of electromagnetic waves with 

some natural surfaces. IRE Trans, on Antennas and Propa. AP-7: 
5324-5329. 

54. Peake, W. H. 1959b. Theory of radar return from terrain. IRE 

Convention Record, 7:27-41. 

55. Poe, G. A., A. Stogryn and A, T. Edgerton. 1971. Determination 

of soil moisture content using microwave radiometry. Final 
Tech. Rep. 1684-1. Aerojet General Corp. El Monte, California. 

56. Priestley, C. H. B. and R. J. Taylor. 1972. On the assessment 

of surface heat flux and evaporation using large-scale para- 
meters. Mon. Weather Rev. 100:81-92. 

57. Richardson, A. J. and C. L. Wiegand. 1977. Distinguishing vege- 
tation from soil background information. Photogram. Eng. and 
Remote Sensing. 48:1541-1552. 

58. Richerson, J. A. 1971, An experimental evaluation of a theore- 
tical model of the microwave emission of a natural surface. 
Tech. Rpt. RSC-27. Remote Sensing Center, Texas A&M University, 
Col lege Station, TX. 

59. Rosenthal, W. D. 1981. Development of visible/infrared/micro- 
wave agricultural classification and biomass models. Ph.D. 
Dissertation. Dept, of Agricultural Engineering, Texas A&M 
University. College Station, TX. 214 p. 

60. Rouse, J. W. Jr. 1972. The effect of the subsurface on the 
depolarization of rough-surface backscatter. Radio Sci. 7: 
889-895. 

61. Rouse, J. W. Jr., H. C. MacDonald and W. P. Waite, 1969. Geo- 
science applications of radar sensors. IEEE Trans. Geosci. 
Elect. GE-7(1) :1-19. 

62. Rouse, J. W. Jr., R. H. Haas, J. A. Schell and D. W. Deering. 

1973. Monitoring vegetation systems in the great plains with 
ERTS. Third ERTS Symp. NASA SP-351. 1:309-317.' 

63. Schmugge, T. J. 1976. Preliminary results from the March 1975 
soil moisture flight. Goddard Space Flight Center Rep. 913-76- 
216. Greenbelt, MD 23 p. 

64. Schmugge, T. J. 1978. Remote sensing of surface soil moisture. 
J. Appl, Meteorology 17:1549-1557. 


103 



65. Schmugge, T. J. 1980a. Effect of texture on microwave emission 
from soils. IEEE Trans, on Geosci. and Remote Sensing GE-18 : 353- 
361. 

66. Schmugge, T. J. 1980b. Microwave approaches in hydrology. 
Photogram. Eng. and Remote Sensing ^6 (4) :495-507. 

67. Schmugge, T, J. 1980c. Soil moisture sensing with microwave 
radiometers. Reprints, 1980 Machine Processing of Remote Sensed 
Data Symp., June, p. 346-354. 

68. Schmugge, T. , P. Gloesen, T. Wilheit, and F. Geiger. 1974. 
Remote sensing of soil moisture with microwave radiometers. J. 
Geophys. Res. 79(2):317-323. 

69. Schmugge, T. J., J. M. Meneely, A. Rango, and R. Neff. 1977. 
Satellite microwave observations of soil moisture variations. 
Water Resources Bull. 13:265-281. 

70. Schwarz, D. E. and F. Caspall. 1968. Hie use of radar in the 
discrimination and identification of agricultural land use. Pro- 
ceedings, Fifth Symp. on Remote Sensing of the Environment. 
Report No. 4864-18-X, Willow Run Laboratories of the Institute of 
Science and Technology The University of Michigan, Ann Arbor, p. 

71. Semenov, B, 1966. An approximate calculation of scattering of 

electromagnetic waves from a rough surface. Radioteknika Elec- 
tronika (USSR) 11:1351-1361. 

72. Sibley, R. 1973. Microwave emission and scattering from vege- 
tated terrain. Tech. Rep. RSC-44 Remote Sensing Center, Texas 
ASM University, College Station, TX. 

73. Simonett, D. S. 1970. Remote sensing with imaging radar: a 
review. Geoforum, No. 2:61-74. 

74. Stogryn, A. 1970. The brightness temperature of a vertically 

structured medium. Radio Sci. 5(12) : 1397-1406. 

75. Taylor, R. C. 1959. Terrain return measurements at X, Ku and Ka 
band. IRE Convention Record, 7:19-26. 

76. Taylor, S. A. and G. L. Ashcroft. 1972. Physical edaphology-the 
physics of irrigated and non-i rri gated soils. Freeman and Co., 
San Francisco. 533 p. 

77. Theis, S. W. 1979. Surface soil moisture estimation with the 
electrically scanning microwave radiometer (ESMR). M.S. Thesis. 
Dept, of Meteorology, Texas A&M University, College Station, TX. 


104 


r 


78. Tsang, L., E. Njoku and J. A. Kong. 1975. Microwave thermal 

emission from a stratified medium with non-uniform temperature 
distribution. J. Appl. Physics 46(2) :5X27-5133. 

79. Ulaby, F. T. 1974. Radar measurement of soil moisture content. 
IEEE Trans. Antenna Propa. AP-22(2)257-265. 

80. Ulaby, F. T., M. C. Dobson, and G. A. Bradley. 1977. If you 

want to remotely sense soil moisture, use a C-band radar. Proc. 
Microwave Remote Sensing Symp. Dec. 6-7, 1977, NASA JSC, Houston, 
TX. pp 51-74. 

81. Ulaby, F. T., G. A. Bradley and M. C. Dobson. 1979. Microwave 

backscatter dependence on surface roughness, soil moisture, and 
soil texture: Part II, Vegetation--covered soil. IEEE Trans. 
Geosci. Elect. GE-17(2) : 33-40. 

82. Ulaby, F. T., G. A. Bradley and M. C. Dobson. 1978. Microwave 

backscatter dependence on surface roughness, soil moisture, and 
soil texture: Part I, bare soil. IEEE Trans. Geosci. Elect. GE- 
16(4) : 286-295 . 

83. Ulaby, F. T. , and P. P. Batlivala. 1976. Optimum radar para- 
meters for mapping soil moisture. IEEE Trans. Geosci. Elect, 

GE-14(2) : 81-93. 

84. Waite, W. P. and K. R. Cook. 1974. Broad spectrum microwave 
systems for remotely measuring ecological parameters. Final 
Rpt. NSF-31515. Department of Electrical Engineering, University 
of Arkansas, Fayetteville, AK. 

85. Wang, J. R., 0. Shi ue, E. Engman, J. McMurtrey, P. Lawlass, T. 
Schmugge, T. Jackson, W. Gould, J. Fuchs, C. Calhoon, T. 
Carnahan, E. Hirschmann, and W. Glazar. 1980. Remote measure- 
ments of soil moisture by microwave radiometers at BARC test 
site. Technical Memorandum 80720, NASA Goddard Space Flight 
Center, Greenbelt, MD. 


APPENDIX A 


if 


DATA PROCESSING PROCEDURE 


106 



DATA PROCESSING PROCEDURES 


The processing of the Guymon and Dal hart data sets was 
accomplished through joint effort at Texas A AM University. Tills 
appendix summarizes the work and can also be found in Blanchard 
(I960), McFarland et aTL (1981), and Rosenthal (1981). 

Scatterometer Processing 

Scatterometer data were collected aboard the NASA C-130 in analog 
form on a 14-track tape. Copies of the tape were later sent to Texas 
ASM Uni vers ity/Remote Sensing Center for Individual processing, which 
consisted of three phases. The initial processing converted the 
analog data to digital and copied the digital data onto 9-track 
magnetic tapes. The second phase processed the digital data using 
software which calculated sigma (a 0 ) values for each look angle at 
given time intervals. Data were processed so that a cell size roughly 
had a length of 25 m for 13.3 GHz, 38 in for 4.75 GHz, 50 m for 1.6 
GHz, and 75 m for 0.4 GHz. The processing software was described by 
Claassen et aK (1979) and Clark and Newton (1979). 

After processing scatterometer data, field start and stop times 
were determined for each frequency and polarization from line plots of 
o° vers.j time, and aerial photographs. Times were adjusted by shift- 
ing the start/stop times at least 0.5 seconds to Insure full scattero- 
meter coverage within the field. The final start and stop times 
defined the field boundary and were used in determining field averages 
for each frequency, polarization, and look angle. Time frames during 
excessive aircraft roll and drift (roll greater than 3.5°; drift 
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greater than 0°) wore noted and data from affected look angles were 
deleted from further analysis. 

No known technique or mechanism was available to calibrate all of 
the scatterometers. Consequently, any temporal variation In o° could 
Indicate either soil moisture, roughness, or vegetation changes, or 
unstable sensors. 

Due to excessive aircraft roll and drift, several look angles had 
to be eliminated at Dal hart and Guymon due to the uncertainty of the 
cell being within the field. At Dalhart, data from only one field had 
to be eliminated— field 16 on 8/18/80. Also, data at 40° and 45° look 
angles off nadir from several other fields on 8/18/80 were eliminated 
due to excessive drift (Table A!). At Guymon, flying conditions were 
much worse; consequently, data from more fields needed to be deleted. 
A complete list of omitted look angles is given in Table A2. Data 
from 8/11, 8/14, and 8/17/78 were most questionable. 

Passive Microwave Processing 

The raw analog data collected aboard the aircraft were converted 
to digital uncorrected brightness temperatures at NASA/Goddard Space 
Fight Center (GSFC), Corrected brightness temperatures were 
calculated from an equation developed at NASA/JSC (O'Neill, 1981): 

, , r 2 (T )(L) 

t b ■ T [T u «ppr> - 1 -V - TiM) - • T„] <A1) 

where t is the transmittance of the radome, e is the emissivity of the 
radome, T u is the uncorrected brightness temperature based on raw 
digital counts, L is antenna cable loss factor, Tl Is an antenna 
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TABLE Al. QUESTIONABLE SCATTEROMETER DATA FOR DALHART 
AS OF 5-6-81 


Date 

Field // 

Questionable Analysis 

8/14/80 

Al) data Is good 


8/16/80 

All data Is good 


8/18/80 

L12 R2 20,8,18 

L12 R2 14 

45° (drift 9°) 

40, 45° (drift 11°) 
All Angles 


Lll R3 16 
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TABLE A2. QUESTIONABLE SCATTEROMETER DATA FOR GUYMON 


Date 



Field # 

Questionable Analysis 

8/2/78 

LI 

R1 

2, 4, 6, 7, 8, 2x, lx 

40°, 45° 

(-8° drift, 2° roll) 


L2 

R1 

10,13, 14, 15,2a, 2x, lx 

45° 

(-9° drift) 


LI 

R2 

2, 4, 6, 7, la, 2x, lx 

45° 

(-9° drift) 


L2 

R2 

15,17,2a 

45° 

(-8° drift) 

8/8/78 

L2 

R1 

17, lx 

all 

angles 


L2 

R2 

2A 

all 

angles 


L4 

R1 

26 

all 

angles 


LI 

R2 

2,6,7 

all 

angles 

8/11/78 

LI 

R1 

6,8,2x 

all 

angles 


L3 

R1 

19,22, lx 

all 

angles 


L2 

R1 

2x, 

all 

angles 


L4 

R1 

24,25,27 

aH 

angles 


LI 

R2 

4,6,7, 1A 

all 

angles 


L3 

R2 

22 

all 

angles 


L2 

R2 

10,17 

45° 

(-4° drift, 4° roll) 




2A, 2X 

all 

angles 


L4 

R2 

24,26,27 

all 

angles 

8/14/78 

LI 

R2 

4 

all 

angles 


L3 

R2 

19 

40°, 45° 

(-8° drift, 3° roll) 


L2 

R2 

13 

45° 

(9° drift) 




10 

40°, 45° 

(9° drift, 3° roll) 


LI 

R3 

all fields 

40°, 45° 

(11° drift) 


L3 

R3 

lx 

all 

angles 


L2 

R3 

13,14 

all 

angles 




15 

4i*° 

(9° drift) 

8/17/78 

L3 

R1 

21,22 

35°, 40° 

,45° (-12° drift) 


L4 

R1 

2x, 24, 25,26, 27 

35°, 40° 

,45° (-12° drift) 


L3 

R2 

21,22 

all 

angles 




lx, 19, 20 

40°, 45° 

(-10° drift) 


L4 

R2 

24, 25,2x 

4 K ° 

(-9° drift) 

8/5/78 

LI 

R1 

2 

40°, 45° 



L4 

R1 

2x 

40°, 45° 



L2 

R2 

2x 

40°, 45° 



L4 

R2 

2x 

40°, 45° 
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These fields were deleted from the MFMR plots due to excessive roll; 
drift was not a factor. 

temperature factor, Tr is the radome temperature factor, r 2 Is an 
internal parameter for each frequency, and T 0 is the self-emission 
of the receiver. The various constants used in the equation were 
determined from flights over homogeneous areas. Once brightness tem- 
peratures were calculated, line plots of Tb versus time were pro- 
duced and field start and stop times were determined from the plots. 
The times defined field boundaries and were used in calculating field 
averages for each frequency and polarization. 

Since the passive microwave radiometer was oriented at a constant 
angle (3° from nadir), any excessive roll would imply questionable 
MFMR data. Consequently, any time the airplane had roll greater than 
3.5°, the field average MFMR data were deleted. Table A3 lists the 
deleted data. With the exception of data from one flight line at 
Guymon (1.6 GHz data on 8/11/78 had highly erratic brightness 
temperatures ), brightness temperatures were quite stable. The highly 
variable brightness temperatures indicated local unmeasured variations 
in the field. Therefore, the following fields at Guymon were deleted 
from further analysis: fields 10, 13, 14, 15 and 17. 

NS001/M 2 S Processing 

Raw data collected on the NASA C- 130 were in analog form. The 
data were converted to digital onto 9-track tapes at NASA/dohnson 
Space Center. Included with the surface data were calibration data 
consisting of digital counts from looks at constant radiance targets 
aboard the plane. The calibration data were then used to convert 
digital counts to radiance. 


Ill 



TA8LF. A3. Guymon and Dal hart Questionable MFMR Data 


Date 

Field // 

* Roll 

8/8/78 

L2 

IU 

IX 

5.3 

8/11/78 

L3 

R1 

IX 

4.9 


LI 

R2 

6 

-5.1 


L4 

R2 

24 

4.9' 

8/14/78 

L2 

R1 

10,17,2a 

5. 4, -8, -5. 6 





respectively 


L4 

R1 

27 

4.9 


L3 

R3 

IX 

-4.8 

8/17/78 

L3 

R2 

22 

5.0 

8/18/78 

LI 

R1 

16 

6.3 
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Since radiance is a function of the solar angle, a correction 
factor was needed before comparing crop radiance differences. All the 
Dalhart data were normalized to August 18, the day with the smallest 
solar zenith angle; Guymon data were adjusted to August 11 zenith 
angle conditions. The correction factor used was 

R i 

r » (A2) 

cos o 

where Ri and R c are the non-normal i zed and normalized radiance 
values, respectively, and o is the solar zenith angle. 

Most of the visible/infrared data were good quality at Dalhart 
and Guymon. One of the exceptions was the excessively noisy water 
absorption bands (channels 6 and 7) on 8/14/80 at Dalhart, Since no 
means was possible to correct the bands, they were eliminated from 
further data analysis. Also, channel 1 data for fields 6, 8, 10, 12 
and 22 were deleted due to unstable calibration. 

With the exception of channel 9 (770-863 nrn) MMS data at Guymon, 
the calibration information proved to be quite stable. Table A4a 
lists the equations used to convert raw digital counts to radiance 
values. Note channel 9 had three different equations applicable at 
different periods of the experiment. 

All of the working NS001 bands had less stable calibration 
information at Dalhart. Table A4b lists the equations used to convert 
digital counts to radiance values. Note several channels had 
different calibration value# on each flight day. 

Calibration of the thermal band proved to be different for Guy- 
mon and Dalhart. The calibration, using the PRT-5 data, showed that 
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ORIGINAL PA®* ss 
OF POOR QUALITY 


TABLE A4. Equations Used to Convert Raw NSQ01/MMS Digital Counts 
(DC) to Radiance Values, R, (watts cm" 2 ster“ 1 ) for 
Guyrnon (a) and Dal hart (b) 


a. channel 4 

7 

8 
9 
9 


b. channel 1 
1 
1 
2 
2 

3 
7 

4 
4 
4 


R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 


10.46x10”^ 
m 


(DC-12) 


9.61xl0" 4 
m 


(DC-13) 


8.14xl0~ 4 

2Td 

6.98xl0~ 4 

6.98xl0” 4 

" njo 

6.98xl0" 4 

TKO 

1.96xl0“ 4 
" 2d7 

1.96xl0" 4 
" T51 

1.96xl0” 4 

70 

4.63xl0“ 4 

210 


•* (DC-14) 

-* (DC-12) (8/2, 8/5, and 8/8) 
(DC-10) (8/11) 

^ (DC- 17) (8/14) 

*{DC-1) (8/14 & 8/16 (Fit 1)) 
^ (DC- 1 ) (8/16 (Fit 2)) 

* (DC-1) (8/18) 

*(DC~21) (8/14 - 8/16) 


4, j^ 10 *(DC-21) (8/18) 

-4 

*(DC-29) (8/14-8/16) 

-4 

. 5j - 61xl0 _- *( DC _ 29 ) ( 8 / 18 ) 

-4 

■ *(DC-9) (8/14-8/6 (Fit 1)) 


11.42x10" 4 4 . 

171 


(DC-9) (8/16 (Fit 2)) 


11.42xl0" 4 * 
— To? 


(DC-8) (8/18) 
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TABLE A4. (Continued) 


OMf-Ai Ir/Vi- «-• 

OF POOR QUALITY 


5 

5 

5 

6 

6 

6 


R * ~i§l Xl ° 4 * (PC-8) (8/14-8/16 (Fit 1)) 
R = 4 -*(DC-9) (8/16 (Fit 2)) 

R = - jg j — ■ -- * (DC-9) (8/18) 

R = - -jll 1 - 4 * (DC- 12 ) (8/16) 

R = 4 * (DC- 12) (8/18) 

R = -- ^ (DC-IG) (8/16 & 8/18) 
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at Guymon the low temperature calibration black box aboard the plane 
was too high while the high temperature calibration black box was 
measuring the proper temperature. This implied that low surface 

temperatures were as much as 5 K too high. At Dalhart, the opposite 
condition occurred. The low temperature calibration box was reading 
the proper temperature while the high temperature calibration box was 
reading 5 K too low, suggesting that high surface temperatures were as 
much as 5 K too low. 

The normalization solar correction factors (cos 6i) for Dalhart 
are as follows: August 14, 5.7; August 16 (flight 1), 2.0; and 

(flight 2), 1.1 and August 18, 1.0. For Guymon, the normal ization 

solar correction factors are August 2, 1.7; August 5, 1.6; August 8, 
5.0; August 11, 1.0; August 14, 1.6 and August 17, 1.6. To normalize 
the two data sets, the Guymon data set required a multiplication 

factor of 1.3 to roughly match the radiance values at Dalhart. 

Soil Moisture 

Each sensor has a different cell size. Consequently, to compare 
data, field averages were determined for each sensor and compared to 
ground data field averages. Unfortunately, in some cases, averaging 
point locations of soil moisture proved not to be a reliable field 
average. For instance, several rows were irrigated and seen by the 

sensors but not sampled within the field. Also rainfall events 
occurred at Guymon between sampling periods--on 8/2 and 8/8/78. An 
attempt was made to correct the soil moisture by adding the amount of 
rainfall or irrigation, assuming complete infiltration. In some 
cases, this correction did a good job. But in the end the question- 


able soil moisture data were deleted from the data set. The fields at 
Guyinon with deleted soil moisture data were for 0/2: 22, 27, 20, 25, 
19, 24, 0/8: lx, 2x, 2, 10 and 8/17: lx, (line 2), 

With the deletions, calibrations, and normalizations the Guymon 
and Oalhart data sets were as complete as possible. 
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ORIGINAL pass » 
OF POOR QUALITY 


DALHART BARE FIELDS 


13,3 W (10 DEO) VS VOL, SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) 


13.3 W <15 DEO) VS VOL. SOIL MOISTURE 



VOL, SOIL MOISTURE <0-2 CM) 


13.3 W (20 DEG) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) 


13.3 W (40 DEO) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) 


0 o o Distil) 
+ + + 015*16 
□ a a 017*16 

* ♦ ♦ 013*20 

• • • 021*22 


LEGENOt FIELD 


FIG. B2 Dalhart 13,3 GHz VV o°vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 
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rtf r;:* 

OF POOR QUAU'iY 


GUYMON BARE FIELDS 

4,73 HH (10 DEO) VS VOL, SOIL MOISTURE 4.75 HU (15 DEC) VS VOL. SOIL MOISTURE 



VOL, SOIL MOISTURE <0-2 CM) VOL. SOIL MOISTURE CO-2 CM) 


4.75 HH (20 DEC) VS VOL, SOIL MOISTURE 4.75 HH (40 DEC) VS VOL. SOIL MOISTURE 



0 5 10 IS 20 25 30 0 5 10 15 20 25 30 


VOL. SOIL HOISTUflE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

V 

LEGEND] FIELD X X X CIO 
V Y V G17 
2 Z 2 02 
AAA 02 X 
« • 4t G2I426 
• • • D641U 

FIG. B3 Guymon 4.75 HH o° vs. volumetric soil^moisture (0-2 cm) at 10, 
15, 20, and 40 degree incident angles. 
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original m 

OF POOR QUALITY 


FIG. 


DALHART BARE FIELDS 


4,70 Mil (tO DEO) VS VOL. SOIL MOISTURE 


4.75 HH (15 DEG) VS VOL SOIL MOISTURE 



VOL, SOIL MOISTURE (0-2 CM) 



VOL. SOIL MOISTURE (0-2 CM) 


4.75 HH (20 DEG) VS VOL. SOIL MOISTURE. 4.75 HH (40 DEG) VS VOL. SOIL MOISTURE 




VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGEND* HELD o> 0 o 013*11 
it + + 015*16 

a □ □ 017*10 

* * * 019*20 

• • • 021*22 

MEG 

B4 Dal hart 4.75 HH cr° vs. volumetric soil moisture (0-2 cm) at 
10, 15, 20, and 40 degree incident angles. 
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GUYMON BARE FIELDS 


4,75 HY (10 DEC) VS VOL. SOIL MOISTURE 4,75 HV (15 DEC) VS VOL, SOIL MOISTURE 



0 S 10 15 20 25 30 0 5 10 15 20 25 30 

VOL, SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CMI 


4,75 HV (20 DEC) VS VOL. SOIL MOISTURE 4,75 HV (40 DEG) VS VOL, SOIL MOISTURE 



0 5 10 15 20 25 30 0 5 10 15 20 25 30 


VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGENOi FIELD X X X G10 
V Y Y 017 
Z Z Z 02 
A A A 02X 
tt « # G2U26 
• • • 064 tu 
REG 

FIG. B5 Guymon 4.75 HV a 0 vs. volumetric soil moisture (0-2 cm) at 
10, 15, 20, and 40 degree incident angles. 
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op . pcgr Q^AU'iiY 

GUYMON BARE FIELDS 

1.8 HH (10 DEG) VS VOL SOIL MOISTURE 1.8 HH (15 DEG) VS VOL SOIL MOISTURE 




VOL, SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 


1.6 HH (20 DEG) VS VOL. SOIL MOISTURE 1.6 HH (40 DEG) VS VOL. SOIL MOISTURE 



0 5 10 IS 20 25 30 • 0 5 10 15 20 25 30 


VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LECENUj FIELD XXXGIO 
v * v 017 
Z Z Z 02 
A A A C2X 
» • • 021*26 
• • • G6«m 

FIG. B7 Guymon 1.6 GHz HH a vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 
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GUYMON BARE FIELDS 

1,6 HV (10 DEG) VS VOL SOIL MOISTURE 1,6 HV (15 DEG) VS VOL, SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 


1.6 HV (20 DEG) VS VOL. SOIL MOISTURE t.6 HV (40 DEG) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGENDi MELD x x X CIO 
Y v Y GI7 
Z Z Z,G2 
AAA G2X 
• « « 021*26 
. • • o C6*1U 

FIG. B8 Guymon 1.6 GHz HV o u vs, volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 
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DALHART BARE’ FIELDS 


1.6 HV (10 DEG) VS VOL SOIL MOISTURE 1.6 HV (15 DEC) VS VOL SOIL MOISTURE 




VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 


1.6 HV (20 DEC) VS VOL. SOIL MOISTURE 


1,6 HV (40 DEG) VS VOL SOIL MOISTURE 




VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGENO: FIELD o © o 013414 
+ + + DI541S 
ODD 017410 

♦ ♦ ♦ D19420 

• • • 021422 

FIG. B9 Dalhart 1.6 GHz HV a 0 vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 
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GUYMON BARE FIELDS 

0.4 HH (10 DEG) VS VOL SOIL MOISTURE 0.4 HH (13 DEG) VS VOL SOIL MOISTURE 



0 5 10 15 20 25 30 0 5 10 15 20 2S 30 


VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 


0.4 HH (20 DEG) VS VOL. SOIL MOISTURE 0.4 HH (40 DEG) VS VOL. SOIL MOISTURE 



0 5 10 15 20 25 30 0 5 10 15 20 25 30 


VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LECENOj FIELD x X X G10 
Y y r G17 
Z Z Z G2 
a a A G2X 

• • « G21426 

• • • G6414 

FIG. BIO Guymon 0.4 GHz HH 0 0 vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 
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DALHART BARE FIELDS 


0.4 HH (10 DEC) VS VOL, SOU MOISTURE 



VOL. SOIL MOISTURE (0-2 CM, 

I HH (20 DEG) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM1 

LEGENOi FIELO 


0.4 HH (15 DEG) VS VOL. SOIL MOISTURE 



' VOL. SOIL MOISTURE (0-2 CM! 
0.4 HH (40 DEG) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE 10-2 CM) ‘ 

ooo 013411 
+ + + 015416 
ODD D17410 

♦ * ♦ 019420 

• • • 021422 


FIG. Bll Dalhart 0.4 GHz HH o° vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 



GUYMON BARE FIELDS 

0.4 HV (10 DEG) VS VOL SOIL MOISTURE 0.4 HV (15 DEG) VS VOL. SOIL MOISTURE 



0.4 HV (20 DEG) VS VOL. SOIL MOISTURE 0.4 HV (40 DEG) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTURE (0-? CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGEND j FIELD x * x G10 
V Y Y GI7 
Z Z 2 G2 
& A A G2X 
• • • G2U26 
... G6*m 

FIG.* B12 Guymon 0.4 GHz HV o° vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 


ORIGINAL PAGC \3 
OF POOR QUALITY 


DALHART BARE FIELDS 


0.4 HV (10 DEC) VS VOL SOIL MOISTURE 0.4 HV (15 DEC) VS VOL SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) . VOL. SOIL MOISTURE (0-2 EH) 


0.4 HV (20 DEG) VS VOL SOIL MOISTURE 0.4 HV (40 DEG) VS VOL SOIL MOISTURE 



VOL. SOIL MOISTURE (0-2 CM) VOL. SOIL MOISTURE (0-2 CM) 

LEGENOj F1EL0 o o © D134HI 
+ + + D15416 
ODD 1)171)8 

* ♦ ♦ 013420 

• • • 021422 

FIG. B13 Dalhart 0.4 GHz HV ct 5 vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 



